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Dynamics of the Globular Cluster System Associated with M49 (NGC 4472): 
Cluster Orbital Properties and the Distribution of Dark Matter^ 

Patrick Cote^, Dean E. McLaughlin^, Judith G. Cohen^ and John P. Blakeslee^ 

ABSTRACT 

Using the Low Resolution Imaging Spectrometer on the Keck I and II telescopes, 
we have measured radial velocities for 196 globular clusters (GCs) around M49 
(NGC 4472), the brightest member of the Virgo Cluster. Combined with published 
data, they bring the total number of GCs with measured radial velocities in this galaxy 
to 263. In terms of sample size, spatial coverage, velocity precision, and the availability 
of metallicity estimates from Washington photometry, this radial velocity database 
resembles that presented recently for M87 (NGC 4486), Virgo's cD galaxy and its 
second-ranked member. 

We extract the projected kinematics of the full sample of GCs and of separate 
subsamples of 158 metal-poor and 105 metal-rich GCs. In agreement with previous 
results for the global GC kinematics based on smaller datasets, we find that the 
GC system as a whole exhibits a slow overall rotation that is due almost entirely 
to a net rotation of the metal-poor GC subsystem alone. In a spatial average, the 
metal-rich GCs shows essentially no rotation. As a function of galactocentric position, 
the metal-poor GCs rotate roughly about the photometric minor axis of M49 and 
at an approximately constant level of ~ 100-150 km s~^ out to i? ~ 2i?efT- The 
metal-rich GC system shows some evidence (at roughly 1-a significance) for weak 
rotation (J7i? ~ 50 km s~^) beyond R ^ 0.5Rcs, also about the galaxy's minor-axis, 
but in the opposite direction from the metal-poor GCs. Outside of i? ~ RcS, the 
line-of-sight velocity dispersion of the metal-poor GCs exceeds that of their metal-rich 
counterparts by ~ 50%. We also note the presence of a well defined grouping of 10 
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metal-rich GCs that are located at opposite poles along the galaxy's major axis and 
which appear to be rotating at nearly 300 km about the minor axis. This grouping 
may be the relic of a past merger or accretion event. 

The dynamics of the GC system are modeled by using published catalogs and 
number counts to define three-dimensional GC density distributions as input to 
a Jeans-equation analysis. We show that the GC radial velocities alone point 
unequivocally — and independently of X-ray observations — to the need for a massive 
dark halo associated with M49 and the Virgo B subcluster around it. We then use 
a mass model for M49/Virgo B, constructed without reference to any GC data and 
described in detail in a forthcoming paper, to infer the orbital properties of the M49 
globulars. The GC system as a whole is shown to be consistent with an almost perfectly 
isotropic velocity ellipsoid. It is more difHcult to draw any firm conclusions on the 
orbital (an)isotropy of the two mctallicity subsamples, owing to the large uncertainties 
in their individual spatial density profiles and to the poorly observationally defined 
kinematics of the metal-rich GCs in particular. 

After M87, M49 is the second elliptical galaxy for which we have been able to 
demonstrate velocity isotropy in the globular cluster system overall, when no division 
based on GC color or metallicity is attempted. Thus, the data for these two galaxies 
lend support the general assumption of isotropy when using GC kinematics to study 
the dark-matter distribution in early-type galaxies. We also compare the kinematic 
properties of the GC system of M49 to those of M87, M31, and the Milky Way: the 
other galaxies for which samples of one hundred or more GC velocities have been 
accumulated. We argue that, contrary to the traditional view of GCs as non-rotating, 
or slowly rotating, systems, rotation may in fact be common byproduct of the formation 
of GC systems. However, the quantitative details of the rotation are still not clear, 
particularly with regard to the question of possible differences between metal-poor and 
metal-rich globulars. 

Subject headings: galaxies: halos — galaxies: clusters — galaxies: individual (M49) — 
galaxies: kinematics and dynamics — galaxies: star clusters 

1. Introduction 

Systems of globular clusters (GCs) can be found in galaxies of virtually all types, from the 
brightest members of rich galaxy clusters to faint dwarfs in loose groups. Recent studies suggest 
that, over a range of roughly ten thousand in total luminosity, about 0.25% of the baryonic mass 
in galaxies takes the form of GCs (Blakeslee, Tonry & Metzger 1997; McLaughlin 1999a; Blakeslee 
1999). Given this ubiquity, and the relative ease with which they can be detected and studied 
in external galaxies, GCs have long been recognized as invaluable probes of galaxy formation 
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and evolution (see, e.g., Harris 1991 and references therein). To date, observational studies of 
extragalactic GC systems have most commonly relied on GC specific frequencies, metallicities and, 
in a few cases, ages, to constrain formation models for GC systems and their host galaxies {e.g., 
Cohen, Blakeslee & Rhyzov 1998; Kundu et al. 1999; Puzia et al. 1999; Gebhardt & Kissler-Patig 
1999; Larsen et al. 2001; Jordan et al. 2002; Cohen, Blakeslee & Cote 2003). However, additional 
insights into the process of galaxy formation may, at least in principle, be gleaned from the 
dynamical properties of GC systems. 

Not only do radial velocity surveys of extragalactic GCs provide a way of tracing out the 
distribution of dark matter within the host galaxy but, with sufficiently large samples, they offer 
the hope of measuring directly the shape of the GC velocity ellipsoid. Not surprisingly, the first 
dynamical studies of an extragalactic GC system focused on the rich GC population of M87=NGC 
4486 [e.g., Huchra k Brodie 1987; Mould, Oke & Ncnicc 1987; Mould et al. 1990). However, 
with radial velocities available for just 44 clusters, only very weak constraints could be placed on 
the mass distribution — and then only if some assumption (i.e., isotropy) were made a priori 
regarding the GC orbital properties (see Merritt &; Tremblay 1993). As Merritt & Tremblay 
show, dauntingly large samples of order one thousand velocities would be required to constrain, 
simultaneously and self-consistently in the absence of any other information, both the velocity 
anisotropy of a GC system and the underlying mass distribution of its parent galaxy. 

Only recently have multi-object spectrographs on 4m- and 8m-class telescopes begun to 
produce anywhere near the large numbers of radial velocities needed to execute such a program 
in any galaxy. Hanes et al. (2001) compiled a radial velocity database of 278 GCs associated 
with M87 after combining new CFHT velocities with published data from Mould et al. (1990), 
Cohen &; Ryzhov (1997) and Cohen (2000). By using an independent mass model for M87 and the 
surrounding "Virgo A" cluster (McLaughlin 1999b), Cote et al. (2001) showed from these data 
that the M87 GC system as a whole has an almost perfectly isotropic velocity dispersion tensor. 
Considered separately, the metal-poor and metal-rich GC subsystems showed some evidence for 
weak tangential and radial biases, respectively. 

Since M87 contains the richest and most thoroughly studied GC system in the Local 
Supercluster, it was the obvious first choice for such studies. However, its unique location — 
at the dynamical center of the Virgo Cluster as traced by both galaxies and intracluster gas — 
complicates the interpretation of the dynamical properties of its GC system. For instance. Cote 
et al. (2001) noted that the apparent rotation of the metal-poor GCs in M87 might be better 
described as a velocity "shear" similar to that exhibited on larger scales by the surrounding 
Virgo cluster galaxies. If correct, this would suggest a connection between the M87 GCs and the 
surrounding population of (mostly dwarf) galax;ies, perhaps due to the slow infall of the material 
onto M87 itself. 

If we are to understand which properties of the M87 GC system are generic to early-type 
galaxies, and which are unique to M87 as a result of its central location in Virgo, then we require 
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djmamical studies of additional GC systems. An appealing second target is that of M49 (= 
NGC 4472). As the brightest member of the Virgo cluster, M49 contains ~ 6000 GCs projected 
within 100 kpc of its center (Lee et al. 1998; McLaughhn 1999a; Rhode k Zepf 2001) — just 
the number expected for an elliptical galaxy with the luminosity of M49 and a "normal" specific 
frequency. M87, by contrast, has an integrated luminosity about 80% that of M49 but contains 
some 13,500 GCs projected within 100 kpc of its center. This factor-of-three difference in GC 
specific frequency between the two galaxies may indicate that some as yet unidentified "second 
parameter" has played a role in the formation and/or evolution of their GC systems (which 
otherwise appear remarkably similar; see Cote 2003). However, the relevance of the observed 
specific frequency difference recently has been called into question by McLaughlin (1999a), 
who showed that the number of GCs associated with M87 and M49 are both consistent with a 
"universal" ratio of total GC mass to total baryonic galaxy mass, ecc = 0.26 it 0.05%, once the 
contribution of X-ray gas mass in each galaxy is included. The principal distinction between the 
two may then be the fact that M87 sits at the gas-rich center of the main Virgo cluster, while M49 
occupies the core of the relatively gas-poor, "Virgo B" subcluster (Binggeli, Tammann &: Sandage 
1987; Binggeli, Popescu & Tammann 1993; Schindler, Binggeli &: Bohringer 1999). 

The first radial velocity measurements for M49 GCs were presented by Mould et al. (1990), 
who reported a global velocity dispersion of 340 ± 50 km s~^ based on a sample of 26 GCs. 
Building on this work. Sharpies et al. (1998) presented radial velocities for an additional 47 GCs, 
found a difference in velocity dispersion between the metal-poor and metal-rich subsystems (320 
and 240 km s~^, respectively), and noted the presence of apparent rotation among the metal-poor 
subsystem. These conclusions were reinforced by Zepf et al. (2000), who measured velocities for 
another 87 GCs and found global velocity dispersions of 356ib25 km s^"*^ and 221ib22 km s~^ 
for the metal-poor and metal-rich subsystems, with projected rotation amplitudes of 101 and 15 
km s~-^ for the two populations. These findings are in stark contrast with those for M87, where 
both GC populations were found to be rapidly rotating, with projected rotation amplitudes of 
roughly 170 km s~^ (Cote et al. 2001). Zepf et al. also calculated the mass distribution around 
M49 implied by the velocity-dispersion profile of the full GC system under the assumption of 
orbital isotropy, and showed this to resemble the gravitating-mass profile inferred from X-ray 
observations (Irwin & Sarazin 1996). 

In this paper, we present 196 new radial velocities for GC candidates in M49 obtained with 
the Keck I and II telescopes. Combining our new radial velocities with those from the literature, 
we isolate a sample of 263 distinct objects whose colors and radial velocities are consistent with 
their being GCs belonging to M49. This database is comparable to that presented recently for 
M87 (Hanes et al. 2001) in terms of sample size, spatial coverage (from about 0.1 to 3 effective 
radii), velocity precision, and the availability of metallicities from Washington photometry. With 
it, we largely confirm the earlier conclusions summarized above on the global, spatially averaged 
kinematics of the GC system as a whole and of its metal-poor and metal-rich components. 
However, with our expanded dataset, we are better able to examine trends in the GC kinematics 
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as a function of galactocentric position. We see some evidence that the bulk of the metal-rich GC 
system may in fact rotate — albeit very slowly and with low statistical significance — in a direction 
opposite to the metal-poor GCs. 

We then go on to carry out a dynamical analysis of the three GC samples (full, metal-poor, 
and metal-rich) based on the Jeans equation following the methodology laid out in McLaughlin 
(1999b) and Cote et al. (2001). We adopt a mass model for M49/Virgo B that we have constructed 
completely independently of any GC data (McLaughlin & Cote 2003), which allows us to treat 
the velocity anisotropy of the GC system as a free parameter to be constrained. Our main new 
result is the demonstration that velocity isotropy suffices fully to describe the velocity ellipsoid of 
the full globular cluster system {i.e., with no attempt made to separate red GCs from blue). M49 
is the second galaxy (after M87) in which we have attempted this sort of analysis, and it is the 
second in which isotropy has been found to hold to a high level of precision. Quite apart from 
any bearing that GC dynamics will have on issues of galaxy and globular cluster formation, this 
suggests that it could very well be reasonable to assume orbital isotropy in other GC systems in 
order to use them as accurate probes of the total mass distribution in early-type galaxies whose 
mass profiles have not been secured through independent means. 

2. Observations, Reductions and the Composite Database 

In this section, we describe the criteria used to target candidate GCs for observation with the 
Low Resolution Imaging Spectrometer (LRIS; Oke et al. 1995), the procedures used to acquire 
and reduce the LRIS spectra, and the methods used to define the composite GC database for M49. 

2.1. Sample Selection 

Before embarking upon a radial velocity survey, we require a sample of objects whose 
magnitudes and colors are consistent with those expected for GCs associated with the program 
galaxy. We used the photometric catalog of Geisler, Lee &; Kim (1996) to select candidate GCs 
for spectroscopic observations. This catalog, which contains GC candidates spread over an area of 
~ 256 arcmin^ centered on M49, has the benefit that the imaging on which it is based was carried 
out in Washington C and Ti filters, a combination that is known to provide excellent mctallicity 
sensitivity for old stellar populations. A total of 1774 GC candidates brighter than Ti = 23 make 
up the catalog. Since accurate coordinates are needed to design the slit-masks, we combined 
the GC pixel positions given in Geisler et al. (1996) with the coordinates reported by Sharpies 
et al. (1998) for 61 GCs to derive astrometric positions for all objects in the Geisler et al. (1996) 
catalog. 
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2.2. Keck Spectroscopy 

Spectra for candidate GCs were obtained in the course of four observing runs at the W.M. 
Keck Observatory during the 1998, 1999 and 2000 observing seasons. The observing log given 
in Table 1 records the identification number, observing dates, grating, central wavelength and 
resolving power for each observing run. A total of 14 LRIS masks were designed, targeting GC 
candidates in the range 1 < (C - Ti) < 2.25 and 20 < Ti < 22.25.^ Individual masks contained 
between 19 and 29 slits, measuring I'.'O in width and having a minimum length of S'.'O. Two 
different instrumental setups were employed. During runs 1 and 4, we used a 600 line mm~^ 
grating blazed at 5000 A and centered at wavelengths of 5200 A and 4860 A, respectively. During 
runs 2 and 3, we used a 600 line mm~^ grating blazed at 7500 A and centered at 8550 A. For 
both configurations, the dispersion was 1.28 A pixel" ^ giving a spectral resolution of 6 A and 
a wavelength coverage of ~2600 A. For each mask, we obtained a pair of exposures ranging 
between 2400s and 3000s, with comparison spectra for Hg-Ne-Ar-Kr-Xe and quartz lamps taken 
before and/or after each exposure. On each night, long-slit spectra for 2-4 lAU radial velocity 
standard stars were obtained during twilight. The extracted, wavelength-calibrated spectra for 
the candidate GCs were cross-correlated against a master template spectrum created from these 
spectra. All reductions were performed within the IRAF environment.^ 

Uncertainties for radial velocities measured from spectra acquired during runs 1-3 were 
derived from the cross-correlation functions using the method of Tonry &; Davis (1979). A 
total of 147 radial velocities were measured during these observing runs; the mean and median 
uncertainties for these measurements are 104 and 97 km s^^. Spectra obtained during run 4 were 
typically of the highest quality, as only the brightest GC candidates were selected during this 
run {i.e., objects with 20.1 < Ti < 21.51). Radial velocities were measured from these spectra 
following the methods outlined in Cohen & Ryzhov (1997), and representative uncertainties of 
25 or 50 km s^^ were assigned based on the repeatability of measurements based on multiple 
strong absorption lines. Including these measurements, the total number of new radial velocities 
presented here is 196, with mean and median uncertainties of 86 and 79 km s~^, respectively. 



2.3. The Composite Database 

Before examining the GC kinematics or carrying out a dynamical analysis, we must combine 
the new velocities with published measurements to isolate a sample of bona fide GCs. In order 
to correct for possible systematic differences between radial velocities measured during different 



®For {C — Ti) — 1.60 — the mean color of our 263 GCs — this magnitude range is roughly equivalent to 
20-5 ^ 22.7 (Geisler 1996). 

'^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under contract to the National Science Foundation. 
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runs, and with different instruments or telescopes, we have transformed the raw velocities from 
April 1998 and April 1999 (runs 1 and 3, respectively) onto the system defined by our April 2000 
observations (run 4) via the relations, 

Vp,i = {1.023)v'p^^ + 114kms-i 

Vp,3 = {OMO)v'p ;i + 46 km s'^ ^ ^ 

which were derived from GCs whose velocities were measured during run 4, and either of runs 1 
and 3. Since there are no objects in common between runs 2 and 4, no correction was applied to 
the former dataset. Likewise, no corrections were applied to the velocities of Sharpies et al. (1998) 
and Zepf et al. (2000), since both datasets show good agreement with our measurements. However, 
we found it necessary to apply a transformation to the velocities of Mould et al. (1990): 

= (0.904)t>; + 215 km . (2) 

Here Vp^M and v'^ denote the corrected and uncorrected velocities from Mould et al. (1990). 
Both the scale and zero-point terms in this relation agree with those found by Hanes et al. (2001) 
from a comparison of Mould et al. (1990) radial velocities for M87 GCs with those measured at 
Keck and CFHT. Figure 1 compares velocities measured during run 4 with previously published 
velocities and those obtained during other Keck runs, after applying the above transformations. 
The dashed line shows the one-to-one relation, while the dotted line indicates the line of best fit. 

After transforming the radial velocities onto a common system, the various datasets were 
merged and a weighted mean velocity, {vp), was calculated from all measurements for each object. 
Unfortunately, given the galaxy's low systemic velocity (vgai = 997 it 7 km s~^ according to the 
NASA Extragalactic Database) and the large line-of-sight velocity dispersion of the GC system 
((Tp ~ 300 km s~^; see §3 below), there is some ambiguity in distinguishing low-velocity GCs from 
foreground Galactic stars. After some experimentation, we decided to discard those objects with 
radial velocities outside of the interval, 

250 < {vp) < 1950 km s'^, (3) 

where the upper limit was chosen to include a few unambiguous GCs with radial velocities of 
{vp) ^ 1900 km s~^. Requiring this range to be symmetric about the systemic velocity would give 
a low-end cutoff of ~ 50 km s~^. Such a cutoff would almost certainly result in the inclusion of 
some halo stars in our sample. If the velocity dispersion of the M49 GC system is ~ 300 km 
then, in a sample of the present size, we expect at most two M49 globulars will fall below the 
low-end velocity limit in equation (3). 

In addition to the above selection on radial velocity, we imposed the condition that true M49 
GCs must have colors in the range, 

1.0 < (C - Ti) < 2.25 mag, (4) 

and adopted a foreground extinction of E{B — V) = 0.022 based on the DIRBE maps of Schlegel, 
Finkbeiner & Davies (1998). According to the relation of Seeker et al. (1995), this extinction is 
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equivalent to E{C — Ti) = 0.045. With the color-metallicity calibration of Geisler 8z Forte (1990), 
this color range translates to a mctallicity interval of —2.15 <^ [Fe/H] <^ +0.8 dex. From optical 
spectroscopy, Beasley et al. (2000) showed the M49 GC system to span a range of —1.6 ^ [Fe/H] 
^ dex, so our selection limits should include the vast majority of GCs in M49. 

Figure 2 presents an illustration of these selection criteria. The upper left panel shows color 
histograms for the 1774 objects with Ti < 23 in the the catalog of Geisler et al. (1996), as well as 
for the 276 GC candidates with measured velocities (upper and lower histograms, respectively). 
The dashed vertical lines show the adopted cutoffs on (C — Ti) color, while the dotted vertical lines 
indicates the value of (C — Ti) = 1.625 used to divide the sample into metal-poor and metal-rich 
components. According to the color-metallicity relation of Geisler &i Forte (1990), this dividing 
point corresponds to a metallicity of [Fe/H] = —0.68 dex.^ 

The upper right panel of Figure 2 shows a radial velocity histogram for the same sample of 276 
objects. Note that two objects (identification numbers #4497 and ^^1982 in Geisler et al. 1996) 
fall outside the plotted region. The vertical lines show the adopted cutoffs on radial velocity, 
while the arrow indicates the velocity of M49 itself. The lower panel of Figure 2 shows the joint 
constraints on color and radial velocity in the form of a color-velocity diagram. Filled circles show 
all objects with measured radial velocities, with the exception of #4497 and #1982 which, once 
again, fall outside the plotted region. The dashed box indicates the adopted selection criteria for 
color and radial velocity. Note that #2256, the lone circled object inside the boxed region, was not 
included in the final sample of GCs since two independent velocity measurements for this object 
differ by more than 700 km s~^. Similarly, two independent measurements for #1982 disagree on 
whether this object is a true GC or a background galaxy. Discarding it, and the 12 objects that 
fall outside of the dashed region in Figure 2, brings the final sample to 263 GCs. 

The final GC database is given in Table 2. The first five columns of this table record the 
object identification number from Geisler et al. (1996), right ascension, declination, distance from 
the center of M49 in arcseconds, R, and position angle, B, in degrees East of North. In calculating 
R and 6, we take the center of M49 to be a(J2000) = 12:29:47.5 and (5(J2000) = +08:00:10.5 
following Sharpies et al. (1998). The next three columns of Table 2 give the Ti magnitude, 
(C — Ti) color, and metallicity based on the color-metallicity relation of Geisler Sz Forte (1990). 
Individual radial velocity measurements and their source are presented in the next two columns. 
The final column gives the weighted mean velocity and uncertainty for each object. The 263 
confirmed GCs are listed first in the table; objects classified as foreground stars or background 
galaxies, and objects with discrepant radial velocity measurements, are given at the end. 

The spatial distribution of the final sample of GCs is shown in Figure 3. Metal-poor GCs are 
indicated by circles, and metal-rich GCs by squares. Open and filled symbols indicate GCs with 



^Globular cluster metallicities quoted in this paper are ~ 0.1 dex lower than those reported in Geisler et al. (1996) 
since we have adopted the Schlegel et al. (1998) extinction of E{B — V) = 0.022; Geisler et al. (1996) assumed 
E{B -V)=0. 
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positive and negative velocity residuals, Aup = (vp) — Vgai, with respect to M49. Symbol sizes are 
proportional to the absolute value of the residual velocity. The dashed lines show the photometric 
major and minor axes of the galaxy, with respective position angles of = 155° and 65° (Kim 
et al. 2000). The large circle shows our determination of galaxy's effective radius, i?efr = 3fl, 
based on our modeling of the Kim et al. surface photometry (see §4.2 below). At our adopted 
distance of 15 Mpc, 1'= 4.363 kpc, so that ReS = 13.5 kpc. The full GC sample spans a range of 
projected galactocentric radius O'A ^ R ^ 9f5, equivalent to 1-9 ^ i? ^ 41 kpc or ~ 0.1-3i?eff- A 
comparison of Figure 3 with Figure 2 of Cote et al. (2001) reveals the azimuthal distribution of 
radial velocity measurements to be somewhat more uniform in M49 than was the case for M87. 

3. Kinematics of the Globular Cluster System 

Before proceeding with a dynamical analysis of the GC system, we examine its kinematics 
in a model-independent way. Our aim is to determine the basic parameters that describe the 
GC system: the average line-of-sight velocity, projected velocity dispersion and, if rotation is 
important, its amplitude and the position angle of the rotation axis. In this section, we consider 
both the global kinematic properties of the GC system and their behavior as a function of 
galactocentric distance. 

We fit the observed line-of-sight velocities of the GCs with the function, 

{vp) = v^ys + {QR) sin(G - Go) , (5) 

where is the projected position angle from a reference axis (taken here to be the North-South 
direction, such that Q is measured in degrees East of North) and Qq is the orientation of the 
rotation axis of the GC system. A complete discussion of the implications of fitting sine curves of 
this type to projected {{vp), @) data is given in Cote et al. (2001). Briefly, in so doing we assume 
that the GC system is spherically symmetric with an intrinsic angular velocity field stratified on 
spheres, and that the GC rotation axis lies exactly in the plane of the sky (in other words, that 
the galaxy is being viewed "edge-on"). In principle, the rotation amplitude, flR, may then be any 
function of galactocentric radius R; fitting equation (5) to our data does not imply an assumption 
of solid-body rotation or of cylindrical symmetry in the velocity field. Our assumption of spherical 
symmetry for the GC system is reasonable given its modest projected ellipticity (0.16 according 
to Lee et al. 1998). 

3.1. Global Kinematic Properties 

Kinematic properties of various subsets of the GC system are summarized in Table 3. The 
first five columns in this table record the radial range spanned by the cluster sample in question, 
the median galactocentric distance, {R), and number, N, of GCs in each sample, the average radial 
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velocity (the biweight "location" of Beers, Flynn, &; Gebhardt 1990), v^, and the rms dispersion 
about this average (the biweight "scale" of Beers et al. 1990), ap. Quoted uncertainties represent 
68% (1-0") confidence intervals, determined from a numerical bootstrap procedure in which 1000 
artificial datasets are individually analyzed after choosing A'^ clusters at random from the actual 
subsample under consideration. 

Rather than let v^ys be a free parameter in the fits of equation (5) to the {vp)-Q data, we 
make use of the expectation (confirmed by the results in Column 4 of Table 3) that the average 
velocity of the GC system should equal that of M49 itself: Vsys = ^^gai = 997 km s~^. In the 
remaining columns of Table 3, we therefore give the position angle of the rotation axis, 6o, the 
rotation amplitude, (Hi?), and the dispersion about the best-fit sine curve, ap^r, obtained with this 
constraint placed on v^ys- 

The first line of each the three subsections of Table 3 presents the results of immediate 
interest: the global kinematics (i.e., those determined using all appropriate GC data with no 
discrimination on the basis of galactocentric position) for our full GC sample and for the separate 
metal-poor and metal-rich subsystems. The remainder of the table applies to further division of 
the GC data into four wide radial bins; these results are discussed in the next subsection. 

The global velocity dispersion of the metal-poor GC system is significantly larger than that 
of the metal-rich GC system: Gp^r = 342 km s"^ versus 265 km s~^. The global velocity dispersion 
of the entire GC system, ap^r = 312 km s~^, is naturally intermediate to these values; it falls 
somewhat closer to the metal-poor value because that sample is ~ 50% larger than the metal-rich 
one. These results are consistent with the previous measurements of Sharpies et al. (1998) and 
Zepf et al. (2000) based on smaller samples. 

Figure 4 plots the individual velocities (vp) against projected position angle Q for GCs at all 
galactocentric radii in each of our three mctallicity samples. A horizontal, broken line in each 
panel of this figure marks the velocity of M49, fgai = 997 km s~^. It is immediately apparent that 
the lower dispersion of the metal-rich GC velocities relative to the metal-poor sample stems in part 
from the absence of even a single GC with (C — Ti) > 1.625 and (vp) > 1500 km — a fact that 
is also reflected in the empty, upper right-hand corner of the box used to identify globular clusters 
in the bottom panel of Figure 2. It would appear that this effect is real: if the metal-rich GC 
velocity dispersion were the same (~ 340 km s~^) as the metal-poor one, then given a sample of 
105 metal-rich clusters with {vp) < 1500 km s~^, we should have also found 8-9 with {vp) > 1500 
km s~^. There is no obvious deficiency in the spatial coverage of the full radial velocity survey 
that might account for the discrepancy (see Figure 3). 

The bold sine curves in the two upper panels of Figure 4 trace the best fits of equation 
(5) to the entire GC dataset and the metal-poor subset when Vsys is held fixed and equal to 
^gai = 997 km s~^ (i.e., with QR and 0o given in Columns 6 and 7 of Table 3). For the full 
sample, the best-fit sine curve has Q,R = 53 km s^^ and &q = 105°. More importantly, we find 
that rii? = 93 km s~^ and 6o = 100° for the metal-poor sample, while is easily consistent with 
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zero for the metal-rich GC system (hence, the fit to it is not drawn on Figure 4). Thus, averaged 

over R <^ 3i?cff in M49, the metal-rich GC system shows essentially no net rotation, while the 
metal-poor GC system does. Our result for the metal-poor rotation amplitude is consistent with 
that of Zepf et al. (2000), who found U.R = 101 km s~^ for a sample of 93 blue GCs. Our finding 
for the metal-rich GC is similarly in keeping with previous analyses (see also Sharpies et al. 1998). 

The position angle of the galaxy's photometric minor axis is Gmm = 65° (Kim et al. 2000) 
and, to within its (large) 1-cr errorbars, our fitted rotation axis for the metal-poor GC system is 
roughly consistent with this. However, it is incompatible at the l.b-a level with alignment along 
the photometric major axis, which is the solution favored by Zepf et al. (2000) for their smaller 
sample of metal-poor GCs (and their full GC system). Though also similar to the photometric 
minor axis of M49, the fitted Oq for the metal-rich GCs is effectively meaningless because of the 
formally null rotation amplitude. 

To summarize these global, spatially averaged impressions, the average velocity of each GC 
sample is in very good agreement with that of M49 itself. The global velocity dispersions of 
the metal-poor and metal-rich samples differ formally by 76 km s~^ (although the 2-cr errorbars 
overlap) with the metal-poor GCs being dynamically hotter than their metal-rich counterparts. 
The GC system appears as a whole to be slowly rotating, but this is purely the result of a net 
signal from the metal-poor GC subsystem alone; the metal-rich subsystem has none. The velocity 
dispersions of the full and the metal-poor GC samples are not changed by the correction for 
rotation, as is apparent from a comparison of the values of ap^r and dp in Table 3. Thus, regardless 
of its statistical significance, rotation is not dynamically important in the M49 GC system. And 
although it is only loosely determined, the average kinematic axis of the metal-poor GC system 
does not differ significantly from the photometric minor axis of M49. 



3.2. Kinematic Properties as a Function of Projected Radius 

Table 3 also presents the GC kinematics that result from further dividing each of the three 
mctallicity samples that wc have defined into four broad radial bins: (1) 25" < R < 150"; (2) 
150" < R< 250"; (3) 250" < R< 350"; and (4) 350" < i? < 570". Although the parameters in 
some bins with small GC numbers must be viewed with caution, the results suggest that there are 
no drastic changes in the average and rms velocities for each sample. 

It is also clear that the lower velocity dispersion of the metal-rich GCs relative to the 
metal-poor GCs, noted above on the basis of globally averaged kinematics, is only apparent 
outside of the first radial bin; the metal-poor and metal-rich clusters within R < 150" ~ O.SReS 
have the same cTp and ap^r- Figures 5 and 6 demonstrate these points visually. Figure 5 is a plot 
of individual GC velocities versus projected galactocentric radii, for the entire sample regardless 
of color and for the metal-poor and metal-rich subsamples separately. Overplotted arc large, open 
squares at the median radius and average velocity of each radial bin, as reported in Columns 2 
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and 4 of Table 3. The vertical errorbars on these large squares represent the dispersion about the 

average velocity (Column 5 of Tabic 3), and the horizontal errorbars delimit the bins themselves. 
The broken horizontal line through each panel shows the systemic velocity of the galaxy. 

The scatter of points about the line {vp) = 997 km s~^ in the top panel of Figure 5 implies 
that the average and dispersion of [vp) in the full GC system are indeed essentially constant as 
functions of galactocentric position — until the region beyond R ^ T.b ~ 33 kpc, where there 
is an obvious dearth of high-velocity GCs in our sample. We suspect that this is an artifact 
of the small sample size, but we cannot rule out with certainty that is a real, physical effect. 
Until the question can be addressed directly (with still more GC velocity measurements at large 
radii), any kinematics referring to i? ^ 30 kpc in the M49 GC system should perhaps be viewed 
as provisional. Meanwhile, the middle and bottom panels of Figure 5 suggest that the average 
velocities of the metal-poor and metal-rich subsamples agree separately, at all i? ^ 30 kpc, with 
the galaxy's systemic velocity. The bottom panel also shows again the absence of any red GCs 
with [vp) > 1500 km s~^, and it further reveals a deficit of clusters with v ^ 700 km s~^ at radii 

We have also constructed "smoothed" , binning-independent radial profiles of average and rms 
GC velocity, following the procedure of Cote et al. (2001): a radial bin of fixed width (chosen here 
to be AR = 120" ~ 8.7 kpc) is slid through the GC dataset, centering on each globular in turn, 
and if the number of objects falling within this spatial bin exceeds a minimum of 20 then the 
biweight average and rms velocity at that position are calculated. (Equation [5] is also fit to the 
{vp)-Q data in the bin around every point; the rotation profiles that result are discussed below.) 
Confidence intervals about the best-fit kinematics are defined by the same bootstrap method that 
was used to calculate the errorbars in Table 3. 

The profile of average GC velocity as a function of galactocentric position that results from 
this procedure only confirms that is constant (inside R 30 kpc) and equal to Vgai for each 
metallicity sample. The velocity dispersion profiles that we obtain are shown in Figure 6 for the 
full GC system and the metal-poor and mctal-rich subsets (filled squares represent the dispersion 
(Tp of velocities about the average; open squares, the dispersion ap^r about the best-fit sine curve 
at each point). The horizontal errorbar in the top panel represents the 2' width of the sliding 
radial bin that we have used; only points separated by this distance in these plots are statistically 
independent. The dotted and solid lines around the points in each panel denote the 68% and 95% 
confidence intervals for each measurement. All in all, the impressions gleaned from Table 3 and 
Figure 5 are confirmed in Figure 6. 

The velocity dispersion of the full GC system is remarkably constant as a function of radius 
inside i? ^ 30 kpc, never straying by more than la from the global average of 312 km s^^ (dashed 
horizontal line in the top panel). We do not see the gradual decline in ap{R) with radius that 
was inferred by Zepf et al. (2000) on the basis of a heavily spatially smoothed representation of a 
smaller dataset. The dispersion profile of the metal-poor GCs is also relatively flat; if anything, it 
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tends to increase towards larger radii (although this trend is not highly significant) . This behavior 
is opposite to that found by Zepf et al. (2000). 

The metal-rich GC system is roughly as dynamically hot as the metal-poor GC sample at 
small galactocentric radii, but it becomes colder beyond i? ~ 3' ~ 13 kpc ~ 1-Reff before possibly 
increasing slightly again towards i? ~ 6' ~ 2i?cff- None of this variation, however, is significant at 
more than the 68% confidence level. Also, note that we are unable to trace the metal-rich GC 
kinematics beyond i? = 6' in the present approach, because there are fewer than 20 red clusters 
between that radius and R ~ 9f5, whereas every point in Figure 6 relies on data for at least 20 
GCs within a 2' interval in R. Thus, we have included in the bottom panel of Figure 6 (as a 
large open square with vertical errorbars corresponding to 68% and 95% confidence intervals) the 
estimate for cTp in the outermost radial bin of Table 3. Zepf et al. (2000) presented a smoothed 
velocity dispersion profile for the metal-rich GC system that appears somewhat flatter than ours 
inside R ^6'. 

The "rotation-corrected" GC velocity dispersions, nowhere, in any of the three 

metallicity samples, dramatically different from the dispersions, (Tp, about the average velocity. 
This is yet another demonstration that rotation is dynamically unimportant in the globular cluster 
system of this galaxy. Nevertheless, it is of some interest to consider the rotation properties of 

the GCs as a function of galactocentric radius, as they appear to be somewhat more subtle than 
suggested by the global view of §3.1, in which the metal-poor GCs show a statistically significant, 
if small, net rotation while the metal-rich ones do not. 

We show in Figure 7 plots of GC velocity versus projected position angle, for each of the four 
radial bins defined in Tabic 3. A broken, horizontal line in each panel again indicates the overall 
velocity of M49, Wgai = 997 km s~^, and the bold sine curves trace the best fits of equation (5) 
to the metal-poor (filled points and solid curves) and metal-rich (open points and broken curves) 
GCs in each bin. (The fit parameters are again those in Columns 6 and 7 of Table 3.) This gives 
some indication of the quality of the sine fits we are able to achieve with the data so finely divided, 
and the extent to which such fits are even feasible given the data. Note that we do not present any 
rotation fit for the metal-rich GCs in the bin 250" < R < 350", for reasons that we shall discuss 
shortly. 

For more detail. Figures 8-10 plot the runs of ftR and Oq with galactocentric radius in 
the full GC system and in the metal-poor and metal-rich subsystems, obtained during the same 
"smoothing" process that resulted in the velocity dispersion profiles in Figure 6. 

Figure 8, which presents the results for our full sample of 263 GCs, is included primarily 
for completeness, as it is somewhat more profitable to consider the metal-poor and metal-rich 
subsamples separately. Worthy of note here, however, are the large uncertainties at every radius 
in the Qq and QR determinations (as in Figure 6, 68% and 95% confidence intervals are indicated 
here), even when our full dataset is analyzed in bins as wide as 2'. Much larger mimbcrs of 
GC velocities will be required if the rotation properties of the system are to be characterized 
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point-by-point in this way with any real confidence. 

Figure 9 shows the rotation axis and amphtude as functions of R in the metal-poor GC system 
only. Inside i? ^ 25 kpc, 6o is essentially constant and generally within la of the photometric 
minor axis of M49 (drawn as the bold, solid lines in the top panel of Figure 9). Indeed, it is made 
clear in this plot that the rotation axis of the metal-poor GC system by and large falls closer to 
the galaxy's minor axis than its major axis. The flR profile is similarly rather fiat, holding roughly 
steady (again, within the l-cr uncertainties) at the level of 100-150 km s~^. It is nonzero almost 
everywhere at better than the 68% confidence level, but less than the 2-a level. The obvious 
exception is at very small radii, R < 2'. The strong negative rotation suggested in the innermost 
radial bin of Table 3, and illustrated in the top panel of Figure 7, is seen here simply not to be 
significant; it is an artifact of the small sample size. 

It could be of considerable interest that the metal-poor rotation axis appears to switch 
suddenly to align with the photometric major axis of M49 (indicated by the bold, dashed lines 
in the top panel of Figure 9) at the largest radii probed by our sample, R ^ 6'. This is clearly 
responsible for a similar, though slightly more muted, trend in the full GC system (Figure 8). 
However, since every point plotted in these graphs incorporates all GC data from within 1' on 
either side of it, the measurements beyond 6' all rely on some velocities from i? > 7' ~ 30 kpc, 
which is the regime that was mentioned as somewhat suspect on the basis of Figure 5. A larger 
dataset is required to confirm this potential major-axis rotation at large radii in the metal-poor 
GC system. 

Figure 10 presents the analysis as applied to the metal-rich GC system. Since the 
"smoothed" profiles cannot be calculated beyond i? ^ 6' in this case, a single large data point 
is used to represent the sine-fit parameters for all metal-rich GCs with 350" < R < 570" 
(18 kpc ^ R 41 kpc) together. At small radii R ^ 4', the axis of GC rotation is poorly 
determined but is consistent, within the 95% confidence bands, with alignment everywhere along 
the photometric minor axis (or coincidence with the formal best fit, Qq = 100°, to the metal-poor 
GC rotation). Notably, and despite the vanishing net rotation of the globally averaged metal-rich 
GC sample, QR at these small radii differs from zero at the ~ 1-a level and is of opposite sign to 
that found for the metal-poor GCs. Thus, there is marginal evidence for some counter-rotation in 
the metal-rich GC subsystem, relative to the metal-poor one, inside R ^ 4'.5 ~ 1.5/2eff in M49. 

The situation appears in Figure 10 to change suddenly around R ^ 4^5, with Qi? becoming 
positive, large, and significant at the > 95% confidence level. The rotation axis also appears to 
come sharply into line with the minor axis of the galaxy. However, the reversion to a negative 
rotation of much lower significance in the outermost bin of our sample suggests that this effect 
may be related to some curiosity in the data between 4' <^ R 6' . Indeed, inspection of the 
metal-rich GC distribution in the third panel of Figure 7 shows immediately that the azimuthal 
sampling of metal-rich GC velocities in this radial range is highly incomplete: all but a few data 
points lie within a couple tens of degrees of the photometric major axis on opposite sides of the 
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galaxy (0 = 155° and 335°). 

The strong "prograde" rotation that suddenly appears in Figure 10 comes entirely from 
a small group of 10 metal-rich clusters, all contained in the thin annulus 302" < R < 337" 
(22.0 kpc < R < 24.5 kpc) and located at opposite poles of the photometric major axis. In {{vp), 
0) space, these GCs trace almost perfectly (with a dispersion of only '^lOO km s~^) a sine curve 
with peak and trough at the position angles of the major axis and an amplitude QR = +300 
km s~^. Because all 10 objects are so closely spaced around their median radius of {R) = 319", 
all 10 are included in every point with 4'.5 < R < in Figure 10, and they are thus entirely 
responsible for the behavior in those portions of the smoothed @o and CIR profiles (recall that 
every point in these profiles is computed using all GC velocities within ±1' of the indicated 
position). Similarly, the same 10 objects single-handedly erase a (weak) net negative-rotation 
signal from the rest of the metal-rich GC system. 

In the top panel of Figure 11, we show the velocity-position curve of the GCs in question. The 
sine fit indicated has Vsys = 997 km s~^ and Oq = 65° fixed, and a fitted amplitude of Oi? = 296 
km {2-a lower limit: 145 km s~^). The ID's of the ten objects from Table 2 are: 677, 830, 
929, 1508, 6905, 7043, 7894, 8164, 8740, and 8890. Aside from their galactocentric radius, there is 
no obvious, common characteristic to distinguish them from any other data, although we do note 
that five of the 10 clusters have super-solar metallicities, and that the one [#6905] falling farthest 
from the fitted sine curve is also the most metal-poor. 

The bottom panel of Figure 11 is an attempt to delineate the rotation curve of the metal-rich 
GC system with the clusters at 302" < R < 337" isolated from the rest of the sample. On the 
basis of the smoothed profiles at i? < 4' and the single-bin datapoints at i? > 6' in Figure 10, we 
fix the position angle of the rotation axis to coincide with the galaxy minor axis for definiteness 
here {@o = 65°). The same smoothing procedure is applied as before, to fit for fiR in equation 
(5), but only using GC data interior to R < 300". The result, shown as the small, solid squares 
between 2' ^ < 4' in Figure 11 (but incorporating data from 1' ^ i? < 5'), can hardly be called 
a profile any longer, as now only the two endpoints are independent. It does, however, give an 
impression of the (small) impact that a particular choice of bin boundaries has on the derived GC 
rotation amplitude. The open square at i? = 5^3 = 23.2 kpc represents the clusters in the top 
panel of Figure 11, while the open square aX R = 7'.25 is at the rotation amplitude fit to the GCs 
with 350" <R< 570" with v^ys = 997 km s"^ and Go = 65° both held fixed. 

As we suggested just above, the removal of the small group of apparently rapidly rotating 
clusters around R = 22-24.5 kpc leaves a metal-rich GC system that, beyond i? ~ 2' 2± 0.6i?eff 

appears to be rotating — albeit slowly, and with less than impressive statistical significance — 
roughly around the minor axis in a direction counter to the metal-poor GC system. Whether this 
feature is real or not — and, indeed, whether the highly contaminating influence of the 10 isolated 
GCs is an unlucky byproduct of sampling statistics or an indication of a real, physical grouping 
of newly recognized objects — can only be decided by the acquisition of further velocity data. 



-16- 



Additional observations of these clusters would certainly be useful, as they may be the relic of 

a past merger or accretion event. In any case, it is clear that globally averaged kinematics are 
bound to miss potentially significant and complex radial variations within a GC system. Any 
arguments about GC and galaxy formation and evolution that draw only on globally averaged 
rotation should be regarded in this light. 

Wc have checked that excluding the aberrant metal-rich GCs does not lead to significant 
changes (no more than a few km s~^ for velocities, or degrees in 60) in the derived kinematics 
of the full GC velocity sample, either in Table 3 or Figure 8. It does, of course, alter the global 
rotation of the metal-rich GC sample from that quoted in Table 3, and the revised numbers are 
given in Table 4 below (§5). Their exclusion also leaves too few GCs in the bin 250" < R < 350" 
to derive a reliable rotation solution there, which is why none is given in Table 3. However, 
the spatial sampling of points around an annulus at a given radius should not affect grossly the 
estimation of a velocity dispersion there, and thus the metal-rich ap{R) profile in the bottom 
panel of Figure 6 (which excludes no data) is likely still reasonable enough — although the points 
between 4' ^ R ^ 6' should obviously be regarded with due caution. 

Finally, we plot in Figure 12 the absolute value of the ratio of rotation amplitude to velocity 
dispersion as a function of radius in the full GC system, the metal-poor subsystem, and the 
metal-rich sample. Note that the latter profile corresponds to the rotation solutions as obtained in 
the bottom panel of Figure 11, rather than Figure 10. For the metal-poor subsystem, a constant 
\^lR\/ap ~ 0.3-0.4 could be taken as an adequate description, within the 1-a uncertainties, of the 
metal-poor GC data at all radii (except perhaps around 7' ~ 30 kpc). Thus, as the low statistical 
significance of all the rotation amplitudes has already suggested, rotation is not an important 
source of dynamical support in the M49 GC system. 

The weighted average for the full GC system, from either of Tables 3 or 4, is \QR\/ap ~ 
0.15-0.2. A similar, though perhaps slightly higher, ratio also applies in the metal-poor GC 
system. And, with only the exception of the 10 extreme clusters at i? ~ 23.2 kpc, the metal-rich 
GC system is consistent everywhere with \QR\/ap 0.1-0.2. In the next section, we therefore 
explore the dynamics of the M49 GC system through a Jeans-equation analysis that ignores 
altogether the presence of rotation. Nowhere should this expose us to uncertainty or error at any 
more than a ~10% level. 



We make the simplifying assumption of spherical symmetry to proceed with a dynamical 
analysis based on the Jeans equation in the absence of rotation: 



4. Dynamical Models 



2/3ci(r) 



'nci{r)crj{r) = -nci(r) 




(6) 



r 
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Here nci(r) is the three dimensional density profile of the GC system; ar{r) is its intrinsic velocity 
dispersion in the radial direction; Pci{i") = 1 — crg{r) /a'^{r) is a measure of its velocity anisotropy; 
and Mioi(r) is the total gravitating mass enclosed within a sphere of radius r. 

Our analysis of the M49 GC system parallels that undertaken for M87 by Cote et al. (2001), 
and it differs from the approach normally taken in dynamical studies of cxtragalactic GC systems 
{e.g., Huchra & Brodie 1987; Mould et al. 1990; Cohen k Ryzhov 1997; Sharpies et al. 1998; 
Kissler-Patig & Gebhardt 1998; Zepf et al. 2000), in which the measured line-of-sight velocity 
dispersions are used to infer the deprojected profile (Tr{r) so as to solve equation (6), under 
the assumption of orbital isotropy [Pc\{r) = 0], for the gravitating mass distribution Mtot(r)- 
We opt instead to determine a mass model for M49 and the Virgo B subcluster a priori (and 
independently of any GC data) and then to use this model, along with three-dimensional GC 
density profiles that fit published number counts, to solve equation (6) for the velocity dispersion 
profile ar{r) under a variety of assumptions on the velocity anisotropy /3ci(r). These model profiles 
are then numerically projected for comparison with the observed (Tp{R) profiles obtained in §3.2, 
in order to examine more directly what range of f3c\{r) is compatible with the data. (The observed 
dispersion profiles are still too noisy to be deprojected and used to solve directly for (ic\{f) with 
any confidence.) Romanowsky & Kochanek (2001) present an analysis of the M87 GC system that 
employs different methodology for the same basic goal. 

For later reference, the relevant projection integrals are 

where the surface density profile Nf.i{R) is related to the three-dimensional nci(r) by 

N,(R) = 2l^n,(r)^^. (8) 

Throughout this Section, we discuss all distances, densities, masses, and luminosities in 
physical units that assume a distance of 15 Mpc to M49. 



4.1. Density Profiles for the GC System 

Figure 13 shows the projected number-density profiles, Nci{R), for GCs of all colors 
(mctallicities), taken from the comprehensive and independent studies of McLaughlin (1999a) 
and Rhode &: Zepf (2001). These authors have already corrected the GC densities to remove 
contamination by foreground stars and background galaxies, and to account for GCs not directly 
counted because they are fainter than the limiting magnitudes of their surveys (this latter 
correction assumes the GC luminosity function, or number of clusters per unit magnitude, to have 
a Gaussian shape with a peak at V = 23.75 mag and a dispersion a = 1.3 mag). Thus, the data 
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points shown here genuinely reflect the total number of GCs per unit area on the sky as a function 
of projected radius in M49. 

For use in equation (6), we treat both sets of data together as a single number-density profile, 

to which we fit the projections of two simple three-dimensional density functions. Shown as the 
bold, solid curve in Figure 13 is the best-fit projection of the density profile suggested by Navarro, 
Frenk, & White (1997) for dark-matter halos, i.e., rid = 'no{r/b)~^{l + r/h)~'^: 

nci(r) = 0.19 kpc-^ (r/11.7kpc)-\l + r/11.7kpc)~^ , (9) 

with (correlated) uncertainties of ±0.07 kpc^^ in the normalization and ±1.5 kpc in the scale 
radius 

The bold, dashed line in Figure 13 shows the best-fit projection of one of the simple family 
of galaxy models developed by Dehnen (1993): nd(r) = no(r/6)-T(l ± r/6)T-^. We assume 7 = 
and find 

nci(r) = 0.40 kpc-^ (1 + r/16.7kpc)-'' , (10) 

with correlated uncertainties of ±0.16 kpc~^ in the normalization and ±2.3 kpc in the scale radius. 

The reduced values of these two fits are comparable. We prefer the Navarro et al. (NFW) 
form of equation (9), because the total baryonic mass density (stars and X-ray gas) at large radii 
in M49 appears to fall off roughly as (McLaughlin & Cote 2003), and we expect the GC system 
to follow this behavior (McLaughlin 1999a). However, the particular form of the density profile in 
equation (10) was chosen specifically for its contrasting asymptotic behavior at both small and 
large galactocentric radii; with it, we can check directly the extent to which the functional form of 
nci(r) used in the Jeans equation (6) might influence our conclusions on the GC dynamics. 

We wish also to model the dynamics of the metal-poor and metal-rich GC subsystems 
separately. To define individual density profiles for them, we make use of the catalog of Lee 
et al. (1998), which includes positions and Washington (C — Ti) colors for ~ 2000 GC candidates 
around M49. Although Lee et al. (1998) themselves presented separate density profiles for the 
metal-poor and metal-rich GCs, the criteria they used to divide their sample by color differ slightly 
from those adopted here, in §2.3. We therefore use their catalog directly to rc-dcrivc the surface 
density profiles of those objects with 1 < (C - Ti) < 1.625 and 1.625 < (C - Ti) < 2.25. We count 
only the points that are also brighter than Ti < 23, to which magnitude the Lee et al. catalog 
is better than 97% complete. Because the catalog includes GCs only out to i? 2± 8' from the 
center of M49, whereas the full GC system is known to extend to i? ^ 20' (Harris &; Petrie 1978; 
Rhode & Zepf 2001), we are unable to measure directly the level of its contamination by stars and 
background galaxies. The exact choice of background is, however, constrained by the fact that the 
corrected number density profile for all colors 1 < (C — Ti) < 2.25 should be consistent with the 
shape of Nci{R) in Figure 13. The total background density, to the limiting magnitude of Ti = 23, 
estimated in this way is Nj, ~ 2.0 ± 0.2 arcmin"^, or 0.11 ± 0.01 kpc~^. We then assume that the 
contaminating backgrounds to be subtracted from the separate metal-poor and metal-rich GC 
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density profiles are 60% and 40% of the total Ni,. This 3:2 ratio is somewhat arbitrary, but it is 
consistent with that predicted by the IAS Galaxy model (Bahcall & Soneira 1980; Bahcall 1986) 
for foreground star counts to Ti = 23, in our color ranges, in the direction towards M49. 

Having thus obtained corrected, metal-poor and metal-rich GC surface density profiles, we 
fit each with projections of the same two three-dimensional functions that we fit to the full GC 
system. The results are shown in Figure 14. The top panel presents fits of the NFW function, with 
an adopted background of iV^ = 2.2 arcmin"^. For the metal-poor and metal-rich GC samples in 
turn, 

n^^{r) (X (r/20.7kpc)-^(l + r/20.7kpc)-2 

ra^^(r) oc (r/5.54kpc)-^(l -Fr/5.54kpc)-^ . (11) 

The dashed curves running through the open datapoints in the top panel of Figure 14 are the 

projections of these functions. The bolder dashed curve is their sum; it agrees reasonably well 
both with the observed density profile of metal-poor and metal-rich GCs combined in the catalog 
of Lee et al. (1998) (filled circles), and with the NFW fit (cq. [9]; the bold, solid line) to the 
combined data of McLaughhn (1999a) and Rhode & Zcpf (2001) in Figure 13. 

The bottom panel of Figure 14 shows our projected fits of the Dehnen (1993) density model 
to the metal-poor and metal-rich GC systems, 

n^P(r) oc (l + r/20.1kpc)-^ 

n^\r) oc (l + r/9.51kpc)-S (12) 

where the adopted background in this case is Nj, = 1.9 arcmin"^. Again, the sum of the two fits 
is shown as the bold, dashed curve. It is fairly consistent with the observed total GC surface 
densities and with the Dehnen-model fit of equation (10). 

Unlike the case for the full GC system in Figure 13, we do not consider the density profiles 
for the GC metallicity subsample to be particularly well constrained. The limited spatial coverage 
of the Lee et al. catalog leaves much room for erroneous extrapolations to larger radii, and the 
lack of direct background density estimates makes the Nci{R) estimates for the metal-poor and 
metal-rich subsystems too uncertain for our taste even at directly observed radii. It is nevertheless 
necessary to specify separate density profiles for the subsamples if their dynamics are to be 
modeled individually. Note that estimates of the absolute normalizations of n^^(r) and ^^^(r) 
are not required in the following analysis, since these cancel out of equation (6). 

4.2. The Need for an Extended Dark-Matter Halo in M49/ Virgo B 

Before going on to constrain the orbital parameters of the GCs in M49, we make a brief aside 
regarding the evidence for dark matter in M49 and the Virgo B subcluster around it, about which 
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various claims have been made over the years. Mould et al. (1990) argued on the basis of radial 
velocities for 26 GCs that the galaxy does indeed contain a dark halo, but they were unable to 
rule out the possibility of a constant mass-to-light ratio due to their small velocity sample and 
the unknown orbital properties of the GCs. Working with a sample of 144 GC velocities, and 
assuming that GCs have isotropic orbits, Zepf et al. (2000) argued that the mass-to-light ratio at 
~ 30 kpc (2-2.5 Res) is at least five times that at a distance of a few kpc, implying the presence 
of substantial dark matter on large spatial scales. 

Irwin & Sarazin (1996) used ROSAT X-ray observations of the hot, gaseous corona around 
M49 to derive a mass profile for M49 extending to i? ~ 100 kpc ~ 8i?efT, leading them to conclude 
that dark matter must be present on spatial scales of R > ReS- Inside R ~ 10-15 kpc (~ Res), 
however, Irwin & Sarazin found mass-to-light ratios that are perfectly compatible with the 
canonical stellar value in elliptical galaxies. Indeed, Brighenti &: Mathews (1997) later decomposed 
Irwin Sz. Sarazin's (1996) mass distribution into separate stellar and dark-matter components (as 
we also do, in §4.3), and argued that any dark-matter halo must contribute negligibly to the total 
mass inside R ^ 7.5 kpc 2± 0.5i?eff- 

Saglia et al. (1993) observed a flattening of the stellar velocity dispersion profile in M49 at 
a much smaller projected radius, R ~ 2-3 kpc ~ 0.2i?efr, and argued that this alone implied the 
presence of an dynamically dominant, extended dark-matter halo. However, the large scatter 
in independent measurements of the stellar kinematics beyond R ^ 1.5 kpc (see, e.g., Caon, 
Macchetto & Pastoriza 2000) rather weakens the conclusion; and it is certainly possible to fit 
self-consistent galaxy models, containing no dark matter, to the stellar kinematics if no other data 
are considered {e.g., Kronawitter et al. 2000). 

Our improved GC velocity dispersion profile from §3 can be used to look at this issue again. 
We consider the hypothesis that light directly traces mass in M49, so that Mtot(?') in equation (6) 
comes from integrating the stellar luminosity density profile and multiplying by a constant, stellar 
mass-to-light ratio. We fit the surface brightness profile of the galaxy taken from Kim et al. (2000) 
— converted from their Washington Ti photometry to Cousins i?-band according to the relation 
of Geisler (1996) — with the projection of the three-dimensional luminosity density profile. 



The projected best fit of equation (13) to the surface-brightness data of Kim et al. (1998) is 
shown in the left panel of Figure 15. For comparison, we also plot the surface brightness profile 
measured by Caon, Capaccioli, &: D'Onofrio (1994), converted from B-hand to Cousins R. The fit 
has parameters 7 = 0.7, a = 2.82 kpc, and Ltot = 1-40 x IO^^Lr^q. It is from this model that we 
derive a (projected) effective radius of i?eff = 3fl ~ 13.5 kpc for M49. Our measurement of the 
galaxy's effective radius is intermediate to the values of ReS = 3'.8 found by Caon et al. (1994) 




(13) 



This slight modification of the models of Dehnen (1993) allows for a more gradual transition from 
the inner power-law behavior, j ^ r~^, to the asymptotic j ^ at large radii. 
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and that of Kim et al. (2000), who obtained i?efr = 2f0 by fitting a de Vaucouleurs law to the 

same data analyzed here. While the precise value of i?efr itself has no bearing on our conclusions, 
we note that the surface brightness profile of M49 is not well described by a dc Vaucouleurs law, 
which is why Kim et al. (2000) restricted their fit to the region 7" < R < 260". As their Figure 8 
shows, the fitted model falls well below the measured brightness profile for R ~ 260", leading to 
likely underestimates of the galaxy's total luminosity and effective radius. 

The right panel of Figure 15 illustrates the three-dimensional stellar mass density profile, 
Ps(r) = Toj(r), with Tq the i?-band mass-to-light ratio. We take its value from McLaughlin &; 
Cote (2003), who show that the stellar kinematics at projected radii 0.1 < i? < 1 kpc in M49 are 
best fit by a mass model with no dark matter if Tq = 5.9 Mq L~^q and the velocity ellipsoid of 
the stellar distribution is slightly radially biased: (3^ = 0.3. We note that this mass-to-light ratio 
corresponds to a i?-band value of 8.1 Mq L^^q, typical of the cores of elliptical galaxies, and that 
this constant-M/ L mass model is very similar to the self-consistent galaxy model of Kronawitter 
et al. (2000). 

We are examining here the hypothesis that light traces mass, and thus the total mass profile 
to be used in the Jeans equation is just the stellar mass 

Ms(r) = To r Attx"^ j{x)dx = ToLtot 

^0 

where, again, Tq is independent of radius in the galaxy and has a value of 5.9 Mq L~^q when no 
additional dark-matter component is invoked. Taking Mtot{r) = Ms{r), then, we solve equation 
(6) for the intrinsic radial velocity dispersion profile of the stars in M49 by replacing nci(r) with 
Psir) oc j{r) from equation (13), and by substituting Ps{r) = +0.3 (McLaughlin &; Cote 2003) 
for Pcii'^)- The projection integrals of equations (7) and (8) are then applied — again with j(r) 
in place of nci(r) — to predict the stellar line-of-sight velocity dispersion in the absence of an 
extended dark-matter halo around M49. 

This model is shown as a bold, dotted line in Figure 16, where it is compared with the stellar 
kinematic data of Davies & Birkinshaw (1988), Sagha et al. (1993), and Caon et al. (2000). It 

does a respectable job of reproducing the observed velocity dispersions between 0.1 ^ R <^ 8 
kpc.^ Although the success relics in part on the considerable scatter among the different datasets 
at R ^ 1 kpc, the main point is that the stellar kinematics in M49, by themselves, provide no 
conclusive evidence for the existence of any extended and massive dark-matter halo. 

The bold, solid line in Figure 16 traces the line-of-sight velocity dispersion profile (Tp{R) 



®The failure of the model stellar (Jp{R) to describe the data at i? ^ 0.1 kpc reflects a poor fit of equation (13) to 
the observed surface brightness profile at /? ^ 2" in M49; a likely departure of the stellar anisotropy from the average 
radial bias /3s(r) = 0.3; and the influence of the nuclear black hole (M, ~ 5 x 10^ M©; Merritt & Ferrarese 2001) at 
radii R ^30 pc. No attempt was made to model these effects in detail, as the complexity in the inner 100 pc of the 
galaxy has no bearing on the larger-scale dynamics of interest to us here. 



(r/a)i/2 




1 + (r/a)V2 





(7 - 27) + (r/a)V^ 
1 + {r/ay/^ 
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derived from equations (6)-(8) for the full GC system, when the constant-M/L mass model 
of equation (14) is adopted, the GC density profile is modeled by the NFW form for nci(r) 
(eq. [9]), and orbital isotropy [Pcli''') = 0] is assumed. It fails by a wide margin to describe the 
observed dispersion profile of the GCs, which is taken from the top panel of Figure 6. The bold, 
short-dashed line is the expected ap{R) profile under the same main assumptions but taking the 
GC density profile to be the Dehnen-model fit of equation (10). The lighter, long-dashed lines 
revert to the NFW GC density profile but assume strongly radially and tangcntially biased GC 
velocity ellipsoids (upper and lower curves). None of these alternate models fares any better in 
accounting for the observed GC velocity dispersion at ^ 6 kpc ^ 0.4-0.5i?cfr. Evidently, the 
mass-to-light ratio must increase with galactocentric radius if the high velocity dispersion of the 
GCs is to be understood. 

Regardless of any assumptions on its velocity anisotropy or uncertainties in its spatial 
distribution, the globular cluster system alone therefore provides incontrovertible evidence, 
independent of any X-ray data, for the presence of an extended dark-matter halo around M49. 
It is worth emphasizing that the GC data alone are unambiguous on this point for spatial scales 

R ^ ReS- 

4.3. Total Mass Profile of M49/Virgo B 

Having established the presence of a dark-matter halo around M49, we briefly describe a 
model for it that is developed and discussed in detail by McLaughlin & Cote (2003). The model is 
constructed to match the available stellar kinematics and X-ray data in M49/ Virgo B. It makes 
no reference to any aspect of the GC system, and it can thus be used as one of two components 
contributing to Mtotir) in the Jeans equation (the other being the luminous component of §4.2) 
as a constraint on the velocity anisotropy of the GC system. 

McLaughlin & Cote (2003) consider a variety of possible analytic descriptions for the 
dark-matter density profile in M49/Virgo B. Among these, the "universal" halo profile of Navarro 
et al. (1997), pdm = K {r/rs)^^{l + r/rs)~'^ (of. eqs. [9] and [11]), is one that provides a good 
match to all available data. We therefore adopt this form here, to specify 

pr 

Jo 

The mass Mtot(r) interior to any (three-dimensional) radius r is then the sum of Mdm{r) and 
Ms(r) from equation (14). The unknown parameters in the total mass profile are the dark-matter 
density normalization K and scale length r^, and the stellar mass-to-light ratio Tq. McLaughlin 
& Cote fix these by first constraining the total i?-band mass-to-light ratio at a radius r = 1 kpc: 
-^tot('' < lkpc)/L(r < Ikpc) = 5.9Mq L]^^^, thereby guaranteeing that the final model will 
adequately reproduce the observed stellar kinematics in the core of M49 (cf. §4.2). This leaves 
only two independent parameters in the mass model, and these are fit by comparing to the X-ray 



r.J l + r r. 



(15) 
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mass measurements of Irwin & Sarazin (1996) and Schindler et al. (1999).^° 

The best-fit mass model with an NFW halo is drawn as a bold, solid line in Figure 17. The 
lighter, dashed line is the mass profile of the dark-matter halo alone, given by equation (15) with 
K = 1.12 X 10~^Mq pc and = 196 kpc. The light, dotted line is the stellar mass distribution of 
equation (14) with a fitted i?-band stellar mass-to-light ratio Tq = 5.7Mq (and the galaxy 

parameters 7, a, and -Ltot already fixed by the fit to the surface brightness profile in Figure 15). 
Explicitly, then, our final model for M49/ Virgo B is: 

Mtot(r) 

Note that the X-ray masses interior to r ^ 10-15 kpc in M49 are fully consistent with the 
mass of the stellar component alone there — although they clearly do not require the complete 
absence of any dark matter at the center of the galaxy, this component does not begin to dominate 
the total mass budget until r ^ 20 kpc ~ 1.5i?eff- Again, it is the GC velocity data that show the 
clearest evidence on smaller spatial scales for the presence of a massive dark halo. 



= Ms(r) + Mdm(r-) 
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4.4. Velocity Isotropy in the GC System 

Figure 18 plots the stellar and GC velocity dispersion data as Figure 16 did, but now with 
model curves computed using Mtot(r) given by equation (16) rather than Ms(r) alone. The mass 
model was constructed by requiring, in part, that it agree with the stellar kinematics, and the bold, 
dotted curve in Figure 18 verifies that this is the case. It is worth noting that our model — which 
now includes the dark-matter halo that the GC data show must be present — favors the stellar 
data of Davies & Birkinshaw (1988) and Saglia et al. (1993) over those of Caon et al. (2000) (the 
small, open squares that scatter well below the dotted curve in the Figure). The bold, solid line in 
this plot is the velocity dispersion profile predicted for the full GC system if its three-dimensional 
density profile nci(r) in equations (6)-(8) is given by the NFW functional fit of equation (9), and 
if velocity isotropy, /3ci(^) = 0, is assumed. The excellent agreement between this model curve and 
the GC data over the full range 6 ^ i? ^ 30 kpc implies that the assumption of isotropy is, in 
fact. essentinJ.hj correct. 



^"The total mass in X-ray gas within a sphere of radius r = 100 kpc centered on M49 is only ~ 1/15 the mass in 
stars and ~ 1/200 that in dark matter, so the gas mass itself is negligible. As it turns out, though, the X-ray gas 
traces the dark-matter distribution throughout M49/Virgo B (McLaughlin & Cote 2003), so its contribution to the 
overall mass budget is actually included in the normalization of Mdm(r). 
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The bold, dashed line in Figure 18 is the model (Tp{R) for the GC system if its density 
profile is taken from the alternate, Dehnen-model fit of equation (10). Isotropy is still assumed. 
This model also compares favorably with the data, falling within the 68% confidence interval 
around the ap estimate at almost every point. We note again that we prefer on other grounds the 
shallower extrapolation of the NFW density profile for the GC system, over the steeper fall-off 
of the Dehnen profile (see §4.1). However, the present comparison of the models shows that our 
inference of velocity isotropy in the M49 GC system does not rely on this particular choice. 

Figure 19 gives a close-up view of the comparison between the velocity dispersion data and 
a number of models for the full GC system. The top panel shows the same data as Figures 18, 
16, and 6, along with the usual 68% and 95% confidence bands. The bottom panel shows the 
projected aperture dispersion profile, i.e., the velocity dispersion of all objects interior to a given 
projected radius R, rather than just those in a narrower annulus. For a model, this cumulative 
spatial average of ap{R) is defined by 



with i?min the projected galactocentric radius of the innermost datapoint in a velocity sample 
(-Rmin = 1-9 kpc for us). Although the interpretation of the coarsely averaged a^p profile can be 
subtler than that of the differential profile ap(R), it has the significant advantage of reduced noise 
at large R. 

The bold, solid line in the top panel of Figure 19 is the same as that in Figure 18: the 
predicted ap{R) profile for the full GC system in the mass distribution of equation (16), if nci{r) is 
given by the NFW profile of equation (9) and /3ci(^) = 0. The same curve in the bottom panel is 
just the result of averaging this model profile with equation (17). The bold, short-dashed curves in 
the two panels of the Figure are analogous to the solid lines but assume the Dehnen-model density 
profile of equation (10) for the full GC system. The lighter, long-dashed curves assume the NFW 
profile for nci(r) and show predictions for velocity anisotropies ranging from extreme radial bias 
to extreme tangential bias: from top to bottom, /da = I — (Tg / = 0.99, 0.5, —1, and —99. These 
anisotropies are included in the models as spatially constant, so that the curves corresponding to 
them should be viewed as representative of more realistic, spatially varying functions Pclif) which 
average out to the constants that we have specified. 

The close adherence of the isotropic models to the GC data is striking, particularly in the 
case of the aperture dispersion profile. The spatial averaging involved there effectively blurs local, 
apparently noise-driven features in the differential ap{R) data (e.g., the bump around i? ~ 22 
kpc), such that the a^p data never stray substantially from the predictions of isotropy for either 
assumed density profile nci(r). The possibility of even modestly radial or tangential average 
anisotropy {i.e., Pel = 0-5 or —1) can be ruled out at the 95% confidence level or better in the full 
GC system. 




(17) 



Aside from its intrinsic importance for models of GC formation and evolution, the finding of 
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such a high degree of consistency with velocity isotropy in the full GC system could be of practical 
use for mcasTirements of the gravitating mass profiles in other early-type galaxies, where data may 
not be available in the variety and quality required for the type of modeling that we described 
in §4.3. M49 is now the second galaxy (the other being M87) in which we have been able to 
show by direct means that an isotropic velocity ellipsoid is a highly accurate description of the 
GC kinematics. The assumption of isotropy in other GC systems may therefore be much better 
justified than before, making it possible to use the Jeans equation (6) with f3ci{r) = — provided 
that nci(r) is also well known — to derive Mtot{r) from a set of velocity-dispersion measurements. 

Finally, we construct dynamical models for the metal-poor and metal-rich GC subsystems 
separately. This is a somewhat riskier proposition, given the noisier dispersion profiles of these 
smaller datasets and the much poorer characterization of the density profiles n^P(r) and n^^(r) 
(Figure 14 in §4.1). With these caveats in mind, we compare the observed and predicted velocity 
dispersion profiles for the metal-poor and metal-rich GCs in Figures 20 and 21. The various curves 
in these figures have meanings analogous to those in Figure 19. The bold, solid curves assume 
orbital isotropy and adopt the NFW fits for the density profiles of the subsystems (eq. [11]), 
while the bold, short-dashed curves assume isotropy and the Dehnen-model density profiles of 
equation (12). The light, long-dashed curves assume NFW density profiles and spatially constant 
anisotropies (3ci = 0.99, 0.5, —1, and —99. 

The isotropic NFW density-model curve in Figure 20 provides quite an acceptable description 
of the observed (Jp{R) of the metal-poor globulars until radii R ^ 20-25 kpc, where the subsequent 
rise in the data might suggest a move to a slight radial bias in the GC orbits. However, the model 
remains within the 68% confidence bands on the measurements nearly until i? ~ 30 kpc, the limit 
of reliability. The observed (Tap data in the bottom panel adhere more closely to the isotropic curve 
at still larger radii. The model assuming a Dehnen-type density profile for the metal-poor GC 
system gives a stronger suggestion of radial velocity anisotropy,^^ although the isotropic case is 
still barely consistent (within the 95% confidence bands) with the observed profile. Because of the 
increased noise in the measured dispersions and, more importantly, the poorly constrained density 
profile, it is difficult to draw any firm conclusions on the orbital properties of the metal-poor 
GCs, beyond the obvious: (1) the hypothesis of isotropy is not far wrong, particularly at small 
galactocentric distances i? ^ 20 kpc ~ 1.5i?efr; and (2) any possible radial bias present at larger 
projected radii is likely to be modest. 

An assessment of the orbital properties of the metal-rich GCs is necessarily more difficult. Not 
only is the extrapolation of their observed density profile very uncertain, but the difficulties posed 
by the 10 rapidly rotating objects concentrated around R ~ 23.2 kpc, as discussed in §3.2, leave 
little in the way of data that can be compared reliably with any model. The top panel shows our 
smoothed velocity-dispersion profile from Figure 6, with those parts relying on any GC velocities 



Given a fixed, observed velocity dispersion for a tracer population in a fixed gravitational potential, a steeper 
density profile always implies a more radial anisotropy through the Jeans equation. 
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from within 302" < R< 337" excised. The most that can be said is that, at i? ~ i?efr = 13.5 kpc, 
the data are once again consistent with velocity isotropy, regardless of which GC density profile 
n^^(r) is assumed. Beyond this, however, the steep drop in ap{R) appears unsupportable in any 
dynamical model that we have constructed. It is true that the outer-radius data point suggests 
some return to understandable behavior in the metal-rich GC system (i.e., 1-cr consistency with 
isotropy for either n^^(r) profile), but the trend of the cumulative crap(< R) profile (towards 
model curves suggesting a strong tangential bias in the orbital distribution) also suggests that the 
low velocity dispersion may be problematic. 

If the full GC system is perfectly isotropic, and the metal-poor component is fairly closely so, 
then it stands to reason that the metal-rich GC system should be close to isotropic as well. Thus, 
its behavior beyond IRcS needs to be much better understood empirically before any dynamical 
sense can be made of it. We noted in §3 (see, e.g., Figure 5) that the low velocity dispersion of the 
metal-rich GCs is the manifestation of an absence of red objects with {vp) ^ 1500 km from our 
velocity survey, and a corresponding lack of metal-rich GCs with (vp) ^ 700 km outside about 
one effective radius in the galaxy. It would be useful to know whether these are artifacts of our 
particular sample or indicative of something more profound. For instance, a velocity dispersion 
dp <^ 250 km at galactoccntric radii R > R^g is even lower than we would expect for the 
stellar light of M49 (cf. Figure 18). Very closely related to this issue is the density profile of the 
metal-rich GC system, which will have to be defined reliably out to much larger radii than we have 
been able to do with existing databases. Even the Dehnen-model density profile that we explored 
here, with its behavior na —>■ at large r, predicts GC kinematics that seem considerably too 
"hot" compared with the metal-rich data. Either current radial velocity surveys have missed a 
number of high-velocity, metal-rich GCs at large galactoccntric radii — perhaps due to some 
observational bias or simple bad luck — or the density of metal-rich GCs system falls very rapidly 
indeed to zero beyond several tens of kpc. Deep, wide-field imaging in a metallicity-sensitive filter 
combination can directly address the latter possibility, and it will be essential to developing any 
deeper understanding of GC dynamics in M49. 

5. Comparing M49, M87, M31 and the Milky Way 

A direct comparison between the kinematic and dynamic properties of GCs in M49 and M87 
seems in order, given the similarity in the respective datascts and the methods of analysis. In fact, 
the GC radial velocities now amassed for these galaxies constitute the two largest radial velocity 
samples yet assembled for pure Population II tracers in early-type galaxies. 

It is worth bearing in mind that that the GC systems of these two galaxies are similar in 
many respects. For instance, both have bimodal metallicity distributions, with peaks at [Fe/H] 
~ — 1.5 and —0.3 dex. The ratio of the number of metal poor to metal rich GCs is very nearly 
the same in both galaxies. The GC formation efficiencies for the composite GC systems and the 
two subsystems are quite similar in both cases {e.g., McLaughlin 1999a; Cote 2003). Likewise, the 
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measured ages of the metal-poor and metal-rich GCs suggest that they are — on the whole — old 
and roughly coeval (Cohen et al. 1998; Puzia et al. 1999; Beasley et al. 2000; Jordan et al. 2002; 
cf. Kundu et al. 1999 and Lee k Kim 2000). 

To some extent, this similarity extends to the dynamical properties of the two GC systems. 
Most notably, in both M49 and in M87, the full GC system has an almost perfectly isotropic 
velocity dispersion tensor {e.g., compare Figure 19 above with Figure 13 of Cote et al. 2001). 
There are, nevertheless, some clear differences in the overall GC kinematics between the two 
galaxies. For instance, the metal-poor and metal-rich subsystems in M87 share an essentially 
identical projected velocity dispersion, while in M49, the metal-rich GC sample appears to be 
significantly colder than the metal-poor component. 

Table 4 summarizes the global kinematic properties of the two GC systems. Note that the 
results for M87 differ slightly from those presented in Cote et al. (2001) since we have re-computcd 
the various kinematic properties after fixing the systemic velocity of each GC subsamplc at the 
velocity of M87 (fgai = 1277 ± 2 km s~^; van der Marel 1994), identical to the approach taken in 
§3 for M49. In addition, we have re-calculated the global quantities for the full GC system and 
the metal-rich component in M49, after removing the 10 rapidly rotating, metal-rich GCs at 22 
^R^ 24.5 kpc. 

As Table 4 shows, it is in their rotation properties that obvious differences between the two 
systems appear. In M87, both the metal-poor and metal-rich GC subsystems show relatively 
rapid rotation, with flR/ap^r ~ 0.45. In M49, however, rotation is somewhat more modest among 
the metal-poor GCs (QR/ap^r = 0.27), and it is very weak indeed in the metal-rich GC system 
{VtR/op^r of order 0.1 and consistent with 0). The orientation of the rotation axes in the two 
systems is also noteworthy. In M87, the metal-rich GCs are found to rotate, at all radii, about 
the galaxy's minor axis. This is also true of the metal-poor GCs beyond R ~ 2i?eff in that 
galaxy, but inside this radius, the metal-poor GCs in M87 show major-axis rotation (see Table 4, 
and Cote et al. 2001). In M49, the metal-poor GCs are consistent with minor-axis rotation 
inside R <^ 2i?cff cind suggest a possible transition to major-axis rotation outside that radius. 
Interestingly, the position angles for the rotation axes of the stellar light in M87 and M49 are 
found to be 00 = -23± 10° and 59 ±3°, respectively (Davies & Birkinshaw 1988). Thus, to within 
the measurement errors, and over the range where the GC and stellar data overlap, the metal-poor 
GCs in both cases seem to trace the rotation of the underlying galaxy. It would be interesting to 
extend the stellar kinematic measurements to larger radii, to see if this trend continues, i.e., if 
the stellar rotation axis also "flips" by 90° in M87 especially. Finally, at the radii in M87 where 
the metal-poor and metal-rich GCs both rotate about the minor axis, they do so in the same 
direction. By contrast, there is some evidence that the metal-rich GC system of M49 may rotate 



^^At present, the best stellar kinematic data for M87 (Sembach & Tonry 1996) extend only to about l.SJiefi. At 
2Rcs, the surface brightness of the galaxy is fiv ~ 23.6 mag arcsec"^, well within the range of modern long-slit 
spectrographs on large telescopes. 
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counter to the metal-poor one. We caution that this result is of low statistical significance, as we 
have a sample of only ~ 100 red GCs spread over an order-of-magnitude range in galactocentric 
position, but it clearly merits further attention. 

These differences between the GC systems of M49 and M87 become all the more puzzling 
when one considers the results for the GC systems of late-type galaxies. To date, in only two 
spiral galaxies have samples of a hundred or more GC velocities been accumulated: the Milky 
Way and M31. In the former case, fl.R/a = 0.56 ± 0.15 for the GC system as a whole, and 
QR/a = 0.32 ± 0.20 and 1.05 ± 0.28 for the metal-poor and metal-rich subsystems (Cote 1999). 
In M31, the global value for all GCs is ClR/a = 0.88 ± 0.09, with the metal-poor and metal-rich 
subsystems having ^R/a = 0.85 it 0.09 and 1.10 ± 0.16 (Perrett et al. 2002). Based on this very 
small sample, it seems that rotation may be more important for the GC systems of disk galaxies 
than for those of ellipticals. Clearly, however, the four GC systems show considerable diversity in 
their rotational properties. It seems that the lone common feature among the GC systems of these 
four galaxies is the tendency for the metal-poor GCs to form a rotating population, albeit with 
widely varying levels of dynamical importance (i.e., the percentage of total kinetic energy stored 
in rotation varies by more than an order of magnitude, from ~ 3% in M49 to ~ 40% for M31). 
Furthermore, the metal-rich GCs of M87, M31, and the Milky Way also form rotating populations 
with Q.R/a ratios that are equal to or greater than those of their metal-poor counterparts. In this 
regard, the metal-rich GC subsystem in M49 may be something of an anomaly — although such a 
conclusion may be premature given the current data. 

The kinematic properties of the GC systems studied to date undoubtedly contain important 
clues to the formation histories of their parent galaxies, although the full implications at this point 
remain unclear. Nevertheless, it seems that the traditional picture of GC systems as non-rotating, 
or slowly rotating, entities — a view that can be traced to the slow rotation exhibited by halo 
GCs outside of the central few kpc in the Milky Way — may be in need of revision. The net 
GC rotation detected in an ever-increasing number of galaxies, which we speculate may be a 
common property of GC systems, probably reflects the outcome of a complex interplay between 
any number of processes that could have shaped the observed GC velocity fields: tidal torques 
from neighboring galaxies, spin-up during gravitational collapse, angular momentum transfer and 
energy dissipation during gaseous mergers, and the conversion of orbital angular momentum to 
rotation during gas-poor mergers and accretions. The extraction of quantitative predictions for 
GC dynamics from models of these processes will be a challenging task — as will be the complete 
and accurate empirical characterization of kinematics to test those predictions. 



6. Summary 

We have reported radial velocities for 196 GCs associated with M49, the brightest member 
of the Virgo cluster. Combining these velocities with previously published measurements brings 
the total number of GCs with measured radial velocities to 263. This sample is comparable to 
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that assembled recently for M87 (Hanes et al. 2001) in terms of size, spatial coverage, velocity 
precision, and the availability of metallicities from Washington photometry. Using these data, we 
have carried out a kinematical and dynamical analysis of the M49 GC system which mirrors that 
presented recently for M87 by Cote et al. (2001). 

Wc confirm previous findings that, when considered in its entirety, the GC system of 
M49 shows a modest net rotation tliat is due almost entirely to the rotation of the metal-poor 
subsystem. Likewise, we verify that, in a global average, the metal-rich subsystem shows essentially 
no rotation. However, when examined in greater detail, there is some weak evidence that the 
metal-rich GCs in this galaxy may, in fact, be slowly rotating, but in the opposite direction to 
the metal-poor GCs. The overall negligible global rotation of the metal-rich component is traced 
to a small sample of 10 clusters which appear to show very rapid minor-axis rotation, in the 
same direction as the metal-poor GCs. These objects merit further attention, as they may be 
the relic of a past merger or accretion event. Given the rather limited sample of metal-rich GCs 
with measured velocities, the significance of these findings must await additional radial velocity 
measurements, but it is abundantly clear that conclusions formed on the basis of globally averaged 
kinematics are likely to miss potentially interesting and complex variations within a given GC 
system. 

A comparison with the rcsTilts for M87, M31 and the Milky Way suggests that rotation may, 
in fact, be generic to GC systems, a turnabout from the traditional view of GC populations as 
slowly rotating systems. Any trends in GC rotational properties with metallicity seem to be 
quite complex, and will likely require GC radial velocity data for an expanded sample of galaxies 
before they can be better characterized and understood. We note, however, that in M49 and M87, 
there is a hint that the orientation of the rotation axes for the metal-poor GCs trace those of the 
underlying galaxy. 

Previous studies of M49, its GC system, and its X-ray halo have produced some conflicting 
claims on the need for dark matter. We have shown that the GC radial velocities and density 
profiles alone now provide unmistakable evidence for a massive dark halo in M49/ Virgo B. 
This result, which is corroborated by analyses of its X-ray halo but is independent of the X-ray 
observations, holds for any choice of GC orbital properties. 

Wc have presented a mass model for M49/Virgo B that satisfies all existing observational 
constraints on the mass distribution from optical and X-ray surface brightness profiles and stellar 
kinematical data. A Jeans-equation analysis of the GC radial velocities and density profiles using 
this mass model reveals the GC system as a whole to be perfectly well described by an isotropic 
velocity dispersion tensor. This is also likely to be true of the separate metal-poor and metal-rich 
subsystems, but definite conclusions must await the measurement of additional radial velocities 
and improved surface density profiles for the subpopulations. Indeed, the poorly constrained 
density distributions of metal-poor and metal-rich GCs on large scales may now be the principal 
obstacle limiting our understanding their orbital properties. In any case, the demonstration that 
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the composite GC system of M49, like that of M87, has an almost perfectly isotropic velocity 
dispersion tensor lends support to the general assimiption of isotropy when using the GC systems 
of early-type galaxies to derive gravitational masses. 
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Table 1. Observing Log 



Run 



(d/m/y) 



Date 




1 



16-17/04/1998 



600/5000 
600/7500 
600/7500 
600/5000 



5200 
8550 
8550 
4860 



870 
1425 
1425 

870 



3 
4 



2 



24/03/1999 
6/04/1999 
30/04/2000 
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Table 2. Radial Velocities for Confirmed and Candidate M49 Globular Clusters 
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Source^ 


(Vp) 








(circsec) 


(deg) 


(mag) 


(mag) 


(dex) 


(km s-l) 




(km s-l) 












Globulcir Clusters 
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Table 2 — Continued 
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Table 2 — Continued 
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8 
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64 
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1 
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12:29:52.92 


8 


00:52.5 
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62.3 


20.60 


1 


32 
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553± 50 


4 
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12:29:37.49 


8 


00:56.7 
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287.3 


21.28 


1 


38 


-1.25±0.10 


1400±131 


1 


1400±131 


5282 


12:30:03.57 


8 


00:56.9 


242.5 


79.0 


21.90 


1 


95 


0.08±0.13 


1356± 68 
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1356± 68 
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Source^ 
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01:21.0 
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3 
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2860 
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7 


58:18.4 


153.1 


137.1 


20.27 


1.21 




-10± 15 


4 


-10± 15 


4497 


12:30:00.58 


8 


00:07.4 


193.7 


91.0 


22.18 


1.33 




6960±146 


1 


6960±146 


5323 


12:29:48.83 


8 


00:59.3 


52.4 


21.4 


20.33 
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1357± 65 


M 


1357± 65 


8096 


12:29:20.93 


8 


04:26.8 


471.1 


303.0 


21.81 


1.68 




75± 51 


3 


75± 51 


9086 


12:29:27.33 


8 


05:59.0 


459.9 


319.3 


21.80 


1.30 




10± 54 
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9228 
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06:16.2 
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346.5 


21.02 
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Discrepant Obiects^ 










1982 


12:29:41.67 


7:57:07.4 


202.8 


205.5 


20.89 


1.01 


-2.02ib0.ll 


648± 43 
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2558± 41 




















25800±150 


4 




2256 


12:29:47.53 


7:57:31.4 


159.1 


180.1 


21.24 


1.90 


0.08ib0.13 


242± 61 


3 


830± 25 



954± 28 Z 



^Key to Source codes: (M) Mould et al. (1990); (S) Sharpies et al. (1998); (Z) Zcpf et al. (2000); (1-4) as given in Table 1. 
^Omitted from final GC sample. 
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Tablc 3. Kinematics of the M49 Globular Cluster System-^ 



R {R) N Op Go^ Up/ 

(arcsec) (arcsec) (km s~^) (km s""^) (deg) (km s~-^) (km s"-*^) 



Full Sample: 263 Clusters with 1.00 < (C - Ti) < 2.25 
25-570 234 263 973l| 313lf 105l|| 53l| 312lf 



25- 


150 


110 


58 


9io;^3 


328; 


f33 
-32 


-1681^1 


1 1 5+l°3 


335 r>^Q 


150- 


-250 


207 


87 


102511 


306; 


f35 
-21 






302+1 


250- 


-350 


300 


69 


9931^9 


294; 


f'll 

-26 






307;t^ 


350- 


-570 


418 


49 


9251^ 


329; 


f50 
-48 




123l^f 


33lli 



25-570 



Metal-poor Sample: 158 Clusters with 1.00 < (C - Ti) < 1.625 

251 158 95711 3451?^ 100l|J 93t|? 342 



25- 


-150 


107 


26 


869; 


150- 


-250 


209 


53 


1017; 


250- 


-350 


296 


48 


985; 


350- 


-570 


415 


31 


878" 



25- 


-150 


110 


32 


9461^8 


334; 


150- 


-250 


196 


34 


1047l^f 


231; 


250- 


-350 


311 


21 


1008^1^ 


208; 


350- 


-570 


435 


18 


9931^^ 


274; 



-1-33 
18 



■72 Q9R+62 c:n-l-45 1 SO+H^ oko+M: 

% 350;;1 96131 162^10^ 324^^1 

51 oqo+48 70-1-30 900-1-116 QCi9-l-42 

-58 '->"J'->_47 '"^-26 "^"^^-81 ^^"^-53 

+77 QKQ-l-eS c;o-l-46 -. 97-I-I88 Q«7+52 

i_58 -aoy-eg ~^°-50 -"-^'-ua •^'^'-70 



Metal-rich Sample: 105 Clusters with 1.625 < (C- Ti) < 2.25 
25-570 202 105 9991^^ 265lf^ 195l|| I2t\l 2Qbt\l 



+33 

^-62 


153ltl 


QO+141 
^''-136 


329li 


-f-46 
-35 


97t^^ 




19911 


1+44 








'-16 








+57 
^-56 


130t|^ 




2771^1 



^AU uncertainties correspond to 68% (l-o") confidence intervals. 
^Systemic velocity held fixed in fit of sine curve: ^^gys = 997 km s~^. 
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Table 4. Global Kinematic Properties of GCs in M87 and M49^ 



M87 M492 





N 






60^ 




N 






00* 








(km s"-*^) 


(km s"-*-) 


(deg) 






(km 


s-i) 


(km s"-*-) 


(deg) 




All GCs 


278 




171+3" 


'"-14 


4^+0 09 


253 


316 


-1-27 
-8 




1091^7 


n 1^+015 

"■^>^-0.08 


Metal-Poor GCs 


161 


397+11 


186+^^ 


61+^'' 
"^-18 


"■^'-0.11 


158 


342 


-H33 
-18 


93l«? 




27+01" 
^•^'-0.11 


R < 2i?eff ^ 


46 




81+i!« 


_ 14+47 


23+°-^* 














Metal-Rich GCs 


117 


365^^1 


~ 3 

155_37 






95 


265 


+34 
-13 






n in+o-2^ 



^AU uncertainties correspond to 68% confidence intervals. 

^The 10 metal-rich clusters between 302" < R < 337" in M49 have been omitted. 

■'The rotation axis of the galaxy has position angle —23° it 10° inside R <^ R.^ff (Davics & Birkinshaw 1988). 
^The rotation axis of the galaxy has position angle 59° ± 3° (Davies &; Birkinshaw 1988). 
SReff = 96" ~ 7 kpc for M87 (de Vaucouleurs & Nieto 1978). 
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Fig. 1. — Radial velocities measured for M49 globular clusters from spectra obtained in April 2000 (i.e., run #4) plotted 
against independent velocity measurements from other datasets. The dashed line shows the one-to-one relation, while the 
dotted line shows the least-squares line of best fit. 
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Fig. 2. — (Upper Left Panel) (C — Ti) colors for all 1774 g lobular cluster candidates from the photometric survey of Geisler, 
Lee & Kim (1996) are shown as the dashed histogram. The solid histogram shows the color distribution for the 276 globular 
cluster candidates with measured radial velocities. The dashed lines indicate our adopted selection criteria on color: 1.00 < 

(C — Ti) < 2.25 mag, while the dotted line shows the color selection used to isolate metal-rich and metal-poor subsamples: 
(C — Ti) = 1.625. (Upper Right Panel) Radial velocity histogram for the 276 candidate globular clusters. The vertical lines 
show our adopted selection criteria on velocity; 250 < v < 1950 km s~^. (Lower Panel) Color versus radial velocity for our 
initial sample of 276 objects (two objects which fall outside the plotted region — #1982 and #4497 — are indicated by the 
arrows) . The dotted region shows the joint selection criteria on color and radial velocity. The lone circled point in this box refers 
to object 2256, which was discarded from the dynamical analysis due to discrepant radial velocity measurements. Excluding 
this object, the dashed box contains a total of 263 globular clusters. 



-45- 



AR« (kpc) 

-40 -20 20 40 




Aa (arcsec) 



Fig. 3. — Spatial distribution of the 263 globular clusters with measured velocities. Circles and squares indicate metal-poor 
and metal-rich globular clusters, respectively (see text for details). For both subsamples, open and filled symbols indicate 
objects having positive and negative velocity residuals, Avp = (vp) — tigai, respectively, where tigai = 997 km . The size 
of the symbol is proportional to the absolute value of the velocity residual. The center of the galaxy is marked by the cross, 
while the solid circle shows our estimate for the galaxy's effective radius, i?(,ff = 3.'1 ~ 13.5 kpc, based on a fit (see §4.2) of 
the surface photometry of Kim et al. (2000). The diagonal lines show the photometric major and minor axes of the galaxy as 
determined by Kim et al. (2000). Exterior to the large, dotted circle at R = 30 kpc, our sample shows an apparent dearth of 
globular clusters with (vp) ^ 1500 km s~^. Conclusions about the kinematics of the globular cluster system beyond this point 
should be considered provisional. 
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Fig. 4. — (Upper Panel) Radial velocity, (vp), plotted against position angle, 0, for the full sample of 263 globular clusters. 
(Middle Panel) (vp) versus for the sample of 158 metal-poor clusters. (Lower Panel) (vp) versus for the sample of 105 
metal-rich globular clusters. In all three panels, the broken horizontal line indicates the galaxy's velocity, fgai = 997 km s~^; 
in the top two, the best-fit sine functions from Tabic 3 arc shown as bold, solid curves. The position angles of the photometric 
minor axis of M49 (Kim et al. 2000) are indicated by the vertical arrows at = 65° and 245°. 
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Fig. 5. — (Upper Panel) Radial velocity, (vp), plotted against projected galactocentric distance, R, for the full sample of 263 
globular clusters. The broken horizontal line shows the velocity of M49. (Middle Panel) {vp) versus R for the sample of 158 
metal-poor clusters. (Lower Panel) (vp) versus R for the sample of 105 metal-rich globular clusters. In each panel, the large, 
open squares are at the median _R and average velocity, v^, of the globular clusters in the four coarse radial bins of Table 3. 
Horizontal errorbars define the limiting radii of the bins, and vertical errorbars represent the velocity dispersion, ap, in each 
bin. 
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Fig. 6. — (Upper Panel) Biweight estimates of velocity dispersion as a function of galactocentric radius in the full sample 
of 263 M49 globular clusters. Solid squares represent the dispersion, Cp, about the average globular cluster velocity at each 
point; open squares show the dispersion, <Tp,r, about the best-fit sine curve describing the variation of {vp) with at each 

radius. Dotted and solid lines around the points delineate the 68% and 95% confidence intervals on the dispersion estimates. 
The dashed horizontal line shows the global value of the velocity dispersion for the sample, and the horizontal errorbar shows 
the width (2' ~ 8.7 kpc) of the sliding radial bin used to compute this smoothed profile. (Middle Panel) As above, but for 
the 158 metal-poor globular clusters. (Lower Panel) As above, but for the 105 metal-rich clusters. Also shown is a large open 
square representing the single bin containing the 18 metal-rich globular clusters beyond R = 350" ~ 25.5 kpc. 



-49- 







100 



200 



300 



2000 - 




1000 ====2^=fc^ 





2000 
1000 


2000 
1000 


2000 



1000 








+ 



* _ _ • 



150^R<250" 










100 200 
(deg) 



300 



Fig. 7. — Radial velocity versus projected position angle for metal-poor and metal-rich globular clusters (filled and open 
points) in each of the four coarse radial bins whose kinematics are summarized in Table 3. The broken, horizontal lines in each 
panel represent the systemic velocity of M49. The bold solid and dashed curves are the sine curves that best fit the data for the 

metal- jjoor and mctal-ricli globular clusters (with parameters taken from Table 3). Vertical arrows in the second and bottom 
panels mark the position angle of the galaxy's photometric minor axis (0 = 65° and 245°). 
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Fig. 8. — (Upper Panel) Projected position angle of the rotation axis of the globular cluster system of M49 (all metallicities), 
as a function of galactocentric radius. The bold, solid lines at Go = 65° and 245° represent the minor axis of the galaxy light; 
bold, dashed lines at ©o = —25° and 155° , the photometric major axis. 68% and 95% confidence bands for the fitted ©o are 

shown as dotted and solid lines around the points. (Lower Panel) Projected amplitude of rotation as a function of radius for 
the full globular cluster sample. The 68% and 95% confidence bands are indicated as in the upper panel, and the horizontal 
errorbaj: has a width of 2': the size of the sliding radial bin used to construct the profiles. 
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Fig. 10. — Same as Figure 8, but for the metal-rich globular clusters. Also shown, as the large open squares, are the best-fit 
Go and QR, with 1- and 2-a errorbars, in the wide, outermost radial bin of Table 3. 
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Fig. 11. — Upper Panel) Detection of strong rotation from 10 metal-rich globular clusters within in a narrow range of 
galactocentric radius, 22 < R < 24.5 kpc. The sine curve drawn has a zero at the photometric minor axis of the galaxy, i.e., 
©0 = 245° (Kim et al. 2000), and an amplitude of 296 km (with 1-a uncertainties of +50 km and —70 km s~^; and 2-a 
uncertainties of +90 km s^-'^ and —150 km s"'^). (Lower Panel) Projected rotation amplitude as a function of galactocentric 
radius in the metal-rich globular cluster system, for fits of the sine curve in equation (5) with both the systemic velocity held 
fixed (at the velocity of M49 itself) and the position angle of the rotation axis forced to coincide with the photometric minor axis 
of the galaxy. The smoothing process of Figure 10 is used only on those globular clusters with R < 302" = 22 kpc. The large 
point aX R = 5^3 and Q,R = 296 km represents the 10 globular clusters isolated in the upper panel, while that at il = 7.'25 
and Q.R = —59 km represents the sine fit with Veye and ©o both fixed in the outermost radial bin, 350" < R < 570". 
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Fig. 12. — (Upper Panel) Ratio of the projected rotation amplitude to the hne-of-sight velocity dispersion for M49 globular 
clusters, plotted as a function of distance from the galaxy center. Filled points take the ratio relative to the dispersion, <7p, 
about the average globular cluster velocity; open points use the dispersion, o"p,r, about the best fit of equation (5) at every 
radius. The horizontal crrorbar shows the 2' width of the sliding radial bin used to "smooth" the data. (Middle Panel) As 
above, but for the metal-poor globular clusters. (Lower Panel) As above, but for the metal-rich globular clusters. The rotation 
amplitude in question is taken from the lower panel of Figure 11 rather than that in Figure 10. The large open square represents 
the 10 rapidly rotating, metal-rich globular clusters in the upper panel of Figure 11. For these objects, |f2-R|/iTp = 1.2; the 
corresponding rotation-corrected value is \QR\/crp^r = 2.9. 




Fig. 13. — Projected number-density profile for the M49 globular cluster system. Filled circles are data compiled by 
McLaughlin (1999a); open squares come from the independent study of Rhode & Zepf (2001). Both profiles already include all 
scalings and corrections to account for background and foreground contamination, and for photometric incompleteness in the 

relevant globular cluster number counts. Solid line shows the best fit of a projected three-dimensional law of the NFW form, 
as indicated on the graph. The broken curve is the best-fitting projection of one of the family of density models developed by 
Dehnen (1993). (See also eqs. [9] and [10] in the text.) 
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Fig. 14. — Fits to the projected number-density profiles of metal-poor and metal-rich globular clusters in M49 (open circles 
and squares), defined by us from the catalog of Lee et al. (1998), for two assumed functional forms of the three-dimensional 
nci(r) (see eqs. [11] and [12]). See the text for a description of the various curves. 
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Fig. 15. — (Left Panel) il-band surface brightness profile for M49 (open squares) derived from the CTi photometry of Kim 
et al. (2000). The dashed curve is the best fit, in projection, of the model in equation 13. Dots show the B-band surface 
brightness profile of Caon et al. (1994) after shifting by the mean galaxy color of [B — R) = 1.52. The arrow shows our estimate 
for the galaxy's effective radius, B,^g = 3fl. (Right Panel) Three-dimensional stellar mass density profile corresponding to the 
model shown in the previous panel, for a spatially constant stellar mass-to-light ratio of To = 5.7Mq ^^^q in the i?-band (see 
§4.3). 
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Fig. 16. — Evidence from the motions of its globular clusters that M49/Virgo B must contain a massive dark halo. Data points 
at small galactocentric radii are stellar velocity-dispersion measurements from Davies & Birkinshaw (1988), Saglia et al. (1993), 
and Caon et al. (2000). Filled circles show our velocity dispersion profile for the full sample of globular clusters, taken from 
the upper panel of Figure 6. The heavy, dotted curve is the stellar velocity dispersion profile that would be produced in a 
self-consistent galaxy model: i.e., one with no dark-matter halo, a constant ij-band mass-to-light ratio of To = 5.9 in solar 
units and a stellar velocity anisotropy /3s = 0.3. The heavy solid line is the velocity dispersion profile predicted for the full 
globular cluster system in such a model, if the density profile n^i{r) is given by equation (9) and an isotropic velocity dispersion 
is assumed. Lighter, long-dashed curves refer to models incorporating strong radial and tangential velocity anisotropy in the 
globular cluster system. The heavy, dashed line is the prediction for an isotropic globular cluster system if nci(r) is given 
instead by equation (10). 
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Fig. 17. — Adopted mass model for M49 and the surrounding Virgo B subcluster, from McLaughlin & Cote (2003) . Dotted 
and dashed curves show the separate mass profiles for the galaxy and the surrounding dark matter halo; their sum is the total 
gravitating mass, Mtot{r), drawn as the bold, solid curve. See equation (16) of the text for details. Points represent the X-ray 
mass measurements used to fit the model: filled circles are from Irwin & Sarazin (1996), with 2-ct errorbars indicated (and one 
point not used in the fit shown as an open circle), while squares are from Schindler et al. (1999) with approximate, factor-of-two 
errorbars attached. 
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Fig. 18. — Comparison of the observed stellar velocity dispersion profile (symbol types as in Figure 16) with the predicted 
profile (heavy, dotted curve) given the mass model of Figure 17 and equation (16), for a stellar velocity anisotropy Ps = 0.3. 
Velocity dispersion data for the full globular cluster system are shown as filled circles, with 68% and 95% confidence bands 

indicated. The bold solid and dashed curves arc predicted crp(R) profiles for the globular cluster system assuming an isotropic 
velocity ellipsoid and adopting the density-profile fits of equations (9) and (10). 
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Fig. 19. — (Upper Panel) Velocity dispersion profile, ap{R), for the full sample of 263 globular clusters along with bootstrap 
estimates of the 68% and 95% confidence intervals (from Figure 6). The heavy solid curve is the predicted profile for clusters 
that are embedded in the mass distribution given by equation (16), follow the density profile of equation (9), and have an 
isotropic velocity dispersion tensor. Lighter, long-dashcd curves vary the assumption on tlic velocity anisotropy; from top to 
bottom, /9cl(^) = 0.99, 0.5, —1, and —99. The bold, short-dashed curve is a model assuming isotropy but using the fit of 
equation (10) to n(,i(r). The horizontal errorbar in the lower left corner represents the 2' smoothing width used to construct 
the empirical profile. (Lower Panel) Aperture velocity dispersion profile for the complete sample of globular clusters. Smooth 
curves are the models of the upper panel, spatially averaged through equation (17) in the text. 
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Fig. 20. — Similaj: to Figure 19, except for the sample of 158 metal-poor globular clusters. Model curves now use the 
density-profile fits for n^^{r) in equations (11) (bold solid and light, long-dashed curves) and (12) (bold, short-dashed curves). 
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Fig. 21. — Similaj: to Figures 19 and 20, but for a sample of metal-rich globular clusters. Objects with 300" < R < 350" 
(21.8 kpc < R < 25.5 kpc) are excluded from the data in the upper panel, though not from the spatially averaged, aperture 
dispersion profile in the lower panel. The large, open square in the upper panel represents the velocity dispersion (with 1- and 
2-cr errorbars) of all metal-rich globular clusters beyond R = 350" = 25.5 kpc. Model curves are analogous to those in the 
previous Figures but employ the density-profile fits for n^^(r) in equations (11) (bold solid and light, long-dashed curves) and 
(12) (bold, short-dashed curves). 



